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ABSTRACT 

Predicting  the  behavior  of  a  qualitatively  described  system  of  solid 
objects  requires  a  combination  of  geometrical,  temporal,  and  physical 
reasoning.  Methods  based  upon  formulating  and  solving  differential  equations 
are  not  adequate  for  robust  prediction,  since  the  behavior  of  a  system  over 
extended  time  may  be  much  simpler  than  its  behavior  over  local  time.  This 
paper  present  a  first-order  logic,  in  which  one  can  state  simple  physical 
problems  and  derive  their  solution  deductively,  without  recourse  to  solving 
the  differential  equations.  This  logic  is  substantially  more  expressive  and 
powerful  than  any  previous  AI  representational  system  in  this  domain. 
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A  Logical  Framework  for  Solid  Object  Physics 

Ernest  Davis 


1.    Introduction 

The  behavior  of  solid  objects  is  an  important  and  natural  domain  for  common  sense 
reasoning.  Solid  objects  are  involved  in  most  terrestrial  natural  phemomena,  and  in  nearly 
every  man-made  artifact.  They  are  familiar  to  every  human  from  infancy,  and  fairly  well 
understood  by  childhood.  So  natural  and  familiar  is  solid  object  behavior  that  it  was  long 
believed  to  be  the  fundamental  type  of  physical  behavior;  hence,  the  efforts  of  physicists 
from  the  Greek  atomists  through  the  beginning  of  this  century  to  explain  physical  phenomena 
of  all  sorts  in  mechanistic  terms. 

There  are  two  major  difficulties  in  building  an  AI  system  that  reasons  about  solid 
objects.  First,  the  physics  of  solid  objects  can  be  quite  complex.  Geometric,  temporal,  and 
physical  reasoning  are  all  involved.  Many  different  aspects  of  geometry  are  involved;  the 
topology  of  objects,  their  gross  shape  and  dimensions,  and  their  differential  surface 
characteristics.  The  motions  even  of  simple  physical  systems  can  take  bewildering  forms; 
consider  a  spinning  penny  settling  down  to  rest  on  a  table  top.  The  gyroscope,  a  very  simple 
physical  system,  not  only  violates  common  sense  predictions,  but  even  baffles  common  sense 
understanding  when  perceived.  Even  with  all  the  tools  of  mathematical  and  physical 
analysis,  simple  problems  can  be  hard  to  solve;  for  example,  the  exact  two-dimensional 
analysis  of  pushing  a  peg  into  a  hole  in  the  presence  of  friction  is  very  difficult.  (J.  Schwartz, 
personal  communication). 

The  second  difficulty  is  that  a  useful  AI  system  must  degrade  well  under  data 
limitations;  it  must  be  able  to  derive  useful  partial  results  from  incomplete  and  imprecise 
problem  specifications.  An  AI  system  must  be  able  to  produce  a  crude  analysis  even  in  cases 
where  exact  information  as  to  object  shapes  or  parameters  is  unavailable  or  unspecified  (as  in 
the  middle  of  a  design  process)  or  too  complex  for  effective  use.  An  AI  system  must  also  be 
able  to  reason  generically  about  classes  of  varying  problems.  For  example,  if  we  build  a 
chicken  fence,  we  want  to  reason  that  no  chicken,  within  the  range  of  chicken  shapes  and 
sizes,  can  go  through  the  fence;  we  do  not  want  to  repeat  the  calculation  for  each  individual 
chicken.  In  this  kind  of  qualitative  reasoning,  human  common  sense  far  excels  conventional 
computational  techniques.  Our  task  is  to  create  an  AI  system  with  the  same  ability  to  apply 
common  sense  to  these  problems. 

This  paper  presents  a  logical  framework  for  qualitative  reasoning  about  solid  objects. 
We  will  present  the  structure  of  a  first  order  language  L  in  which  the  physics  of  solid  objects 
can  be  described;  we  will  define  the  semantics  of  L  in  terms  of  a  formal  model;  and  we  will 
demonstrate  the  usefulness  of  L  by  showing  that  interesting  problems  can  be  solved 
qualitatively  by  inference  from  plausible  axioms  expressed  in  L.  The  language  L  is  much 
more  expressive  and  supports  much  richer  inferences  than  any  previous  representation 
scheme  in  this  domain. 

In  concentrating  on  the  representation  and  formulation  of  knowledge,  and  postponing 
any  questions  of  algorithms  or  control  structure,  we  are  following  the  methodology  advocated 
by  Hayes  (1978a).  Despite  recent  criticisms  of  this  approach  (.McDermott,  forthcoming),  it 
still  seems  the  most  promising  avenue  for  research  in  knowedge  representation  and  common 
sense  reasoning.  Indeed,  the  success  of  our  theory  confirms  the  usefulness  of  the  method  and 
shows  that  it  is  possible  to  identify  an  interesting  body  of  physical  reasoning  where 
McDermort's  objections  to  a  logicist  approach  do  not  apply. 
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However,  our  theory  falls  short  of  Hayes'  vision  in  several  respects.  First,  it  is  not,  by 
anv  means,  a  naive  theory.  Rather,  we  have  used  the  entire  machinery  of  Newtonian 
mechanics,  including  concepts,  such  as  total  mechanical  energy,  that  are  sophisticated  and 
have  no  common  sense  analogue.  The  reason  to  do  this  is,  simply,  that  it  is  very  difficult  to 
develop  an  alternative  mechanics  that  is  consistent  and  powerful  and  that  does  not  entail 
absurd  consequences.*  Second,  some  of  the  axioms  that  we  shall  use  may  not  be  true,  even  in 
our  restricted  domain;  they  may  not  even  be  consistent.  We  believe,  however,  that  they  are 
true  and  consistent  "for  all  practical  purposes'.  That  is,  in  practice,  an  inference  engine  using 
these  axioms  will  never  (hardly  ever)  discover  their  inconsistency  or  be  led  to  serious  error. + 
Third,  it  is  at  least  doubtful  whether  we  have  satisfied  Hayes'  dictum  that  obvious  facts 
should  have  short  proofs.  It  depends,  of  course,  on  what  axioms  you  take  as  the  starting 
point,  and  what  you  consider  short.  Readers  may  make  their  own  judgement  based  on 
sections  6  and  7. 

It  should  also  be  noted  that  the  mathematics  used  here  is  not  "qualitative"  in  the  narrow 
s  nse  of  much  recent  work  on  naive  physics  (Bobrow.  85);  that  is,  it  does  not  represent 
quantities  purely  in  terms  of  order  relationships  and  constants.  I  do  not  consider  this  a 
weakness,  however;  as  I  have  argued  elsewhere,  (Davis,  forthcoming),  this  is  an 
unacceptably  restricting  limitation  for  many  types  of  common  sense  reasoning.  The  logical 
system  is  qualitative  in  the  sense  that  it  can  work  with  partially  specified  information.  This  is 
the  only  sense  that  matters. 

In  the  remainder  of  the  paper,  we  develop  our  logical  theory,  and  we  use  it  to  solve 
variants  of  one  particular  kind  of  problem:  the  problem  of  predicting  what  happens  when  one 
or  two  dice  are  dropped  inside  a  funnel  (Figure  1).  We  will  find  that  different  forms  of  this 
problem  involve  a  rich,  interconnected  body  of  geometric  and  physical  knowedge  for  their 
solution. 

Section  2  deals  with  previous  AI  work  in  this  domains.  Section  3  presents  several 
variants  of  the  "die  in  the  funnel"  example,  and  discusses  how  they  can  be  solved  in  our 
theory.  Section  4  discusses  the  geometrical,  temporal,  and  physical  theories  which  underly  the 
theory.  Sections  5.  6,  and  7  are  for  true  believers  in  formalism.  Section  5  is  a  formal 
definition  of  the  ontology  of  the  language  L.  Section  6  presents  the  predicates  and  axioms 
needed  to  solve  one  particular  instance  of  the  die  in  the  funnel  example.  Section  7  carries  out 
the  complete  analysis  of  this  example.  Section  8  gives  our  conclusions. 

2.    Background 

Despite  the  interest  and  importance  of  solid  object  physics  for  common  sense  reasoning, 
AI  research  in  this  area  has,  to  date,  been  quite  limited  —  notably  less  extensive  than 
research  in  analysing  electronic  circuits.  For  example.  Fahlman's  BUILD  program  (1974) 
determined  the  stability  of  a  tower  of  polyhedral  blocks.  De  Kleer's  NEWTON  (1975) 
predicted  the  behavior  of  a  point  mass  sliding  on  a  constraint.  Bundy's  MECHO  (1978)  used 
force  analysis  and  conservation  laws  to  make  physical  predictions  in  situations  of  specialized 
format.  Forbus'  FROB  (1979)  predicted  the  behavior  of  a  point  mass  flying  among 
constraints.  Funt's  WHISPER  (1980)  predicted  the  behavior  of  a  collection  of  object  by 
simulating  it  in  a  retina-like  image.  Novak's  ISAAC  (1981)  identified  English-language 
programs  of  fixed  form,  and  applied  special  case  equations  to  them.  Shoham  (1985)  has 
systematized  the  rules  that  govern  the  local  mobility  of  an  object  within  constraints. 

All  these  programs  have  provided  valuable  insights  into  the  nature  of  physical  reasoning 
about  solid  objects.  As  compared  to  the  power  of  human  physical  reasoning,  however,  they 
are  extremely  limited  in  their  geometrical  expressivity  and  in  the  range  of  physics  they 
understand.  Of  these  systems,  only  BUILD  deals  with  three-dimensions;  and  only  MECHO 


"    McCloskey,  83)  is  an  inieresting  analysis  of  the  theories  ol  monon  actually  held  by  naive  subjeas. 
*  The  inconsistenc\'  of  the  axioms  does  not  make  formal  nonsense  of  the  language,  since  the  semanacs 
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deals  with  the  motion  of  extended  objects.  Many  kinds  of  physical  interactions  are  omitted; 
for  instance,  none  of  these  systems  could  solve  the  simple  problem,  "What  will  happen  if  a 
spinning  ball  is  dropped  onto  the  ground?" 

Another  limitation  of  these  programs  is  more  subtle,  but  more  fundamental.  Almost 
without  exception,  these  programs  are  based  entirely  on  extrapolating  differential  behavior. 
In  order  to  make  a  prediction,  the  program  first  determines  how  each  state  of  the  system  will 
tend  to  change,  and  then  extrapolates  these  changes  to  predict  a  continual  trend  of  change  up 
to  point  that  the  structure  of  the  system  changes.  This  extrapolation  may  be  done  either 
qualitatively,  as  in  FROB  and  NEWTON  or  symbolically,  as  in  MECHO,  or  using  point-by- 
point  simulation,  as  in  WHISPER,  or  by  numerical  integration,  as  proposed  by  McDermott 
and  Bernecky  (personal  communication).  The  same  limitation  applies  to  most  work  on  other 
types  of  physical  reasoning,  such  as  in  (de  Kleer  and  Brown,  85),  (Forbus,  85),  and 
(Williams,  85),  but  there  it  may  be  less  pernicious.  The  major  exception  was  Hayes'  work  on 
liquids  (1978b),  which  handled  extended  time  periods  through  the  use  of  histories. 

For  example,  FROB  (Forbus,  79)  predicts  the  behavior  of  a  bouncing  ball  in  a  well  by 
dividing  physical  space  into  4  regions  (the  interior  of  the  well,  the  bottom,  and  the  two 
sides),  and  dividing  the  velocity  space  of  cne  ball  into  nine  (motionless,  up,  down,  left,  right, 
and  the  four  quadrants.)  (Figure  2)  There  are  thus  27  possible  states  of  the  system.  (4x9- 
9  impossible  states).  The  laws  of  physics  are  then  used  to  determine  which  transitions 
between  states  are  allowed,  and  thus  a  transition  graph  of  states  is  developed.  The  final 
prediction  is  that  the  system  follows  a  path  in  this  transition  graph. 

However  extrapolation  is  done,  this  kind  of  simulation  is  inherently  inadequate  for 
robust  prediction  in  this  domain.  In  any  analysis  of  this  kind,  each  different  set  of  the 
boundary  conditions  that  can  apply  in  different  situations  must  be  considered  a  different" 
system  state.  Each  such  state  must  be  separately  detected,  categorized,  and  analyzed,  and  the 
system's  progress  through  these  states  must  be  systematically  recorded.  In  many  problems, 
however,  such  a  categorization  is  very  difficult,  and  entirely  pointless.  Consider  the  problem 
in  figure  1;  a  small  die  is  released  inside  a  large  steep  funnel.  Many  states  are  possible  in  this- 
system:  the  die  may  be  in  free-fall;  it  maj'  be  colliding  or  in  continuous  contact  with  the  top 
or  bottom  part  of  the  funnel  funnel  on  any  of  eight  vertices,  twelve  sides,  or  six  faces;  it  may 
be  spinning,  sliding,  or  rolling,  up,  down,  or  around  the  funnel.  If  we  wish  to  determine 
which  face  will  end  up  on  top,  we  will  have  to  trace  the  progress  of  the  die  through  all  these 
states.  Such  a  tracing  is  known  to  be  difficult.  However,  if  all  we  wish  to  determine  is  that 
the  die  falls  through  the  funnel,  then  all  these  intermediate  state  should  not  matter,  and  it 
should  not  be  necessary  to  trace  them.  To  determine  by  any  kind  of  simulation  that  the  die 
falls  through  the  funnel  will  be  almost  as  difficult  as  to  determine  how  it  lies.  This  seems 
implausible. 

Two  further  arguments  may  be  brought  against  simulation.  First,  in  a  problem  like  the 
die  in  the  funnel,  the  sequence  of  states  traversed,  and  the  final  lie  of  the  die,  depend  very 
delicately  on  the  exact  shapes,  sizes,  and  physical  properties  of  the  die  and  the  funnel.  On  the 
other  hand,  the  conclusion  that  the  die  comes  out  the  bottom  is  quite  robust  with  respect  to 
small  variations  in  these  parameters.  Therefore,  if  the  problem  is  specified  with  some  small 
degree  of  imprecision,  simulation  will  either  be  impossibile,  or  involve  some  monstrously 
branching  tree  of  possibilities.  But  in  qualitative  reasoning  we  would  like  the  prediction  that 
the  die  comes  out  the  bottom  to  be  almost  as  easy  with  imprecise  data  as  with  precise. 
Second,  the  complexity  of  simulation  goes  up  rapidly  with  the  number  of  interacting  objects. 
In  figure  3,  for  example,  with  one  die  inside  another  dropped  inside  a  funnel,  the  set  of 
system  states  is  the  cross  product  of  the  possible  interactions  of  the  two  dice  with  the  possible 
interactions  of  the  outer  die  and  the  funnel.  Nonetheless,  the  prediction  that  the  two  dice 
will  come  out  the  bottom  is  almost  as  easy  as  in  figure  1,  with  only  one  die.  (These 
arguments  are  given  at  greater  length  in  [Davis,  86b].) 

In  short,  formulating  and  solving  differential  equations  is  an  inadequate  technique  in 
this  domain,  since  the  behavior  of  these  physical  systems  over  extended  time  is  often  easier 
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to  characterize  than  their  behavior  over  local  time.  A  powerful  physical  reasoning  system 
must  be  able  to  ignore  these  multitudes  of  trees  to  focus  on  the  forest.  We  must  be  able  to 
infer  the  general  quality  of  a  course  of  events  from  broad  characterizations  of  the  physical 
properties  of  the  objects  involved,  without  calculating  each  little  subevent. 

The  programs  cited  do  use  some  techniques  other  than  simulation.  MECHO  and 
NEWTON  use  energy  conservation  to  prune  possible  system  behaviors.  Any  state  with  more 
mechanical  energy  than  the  starting  state  can  be  ruled  out  as  a  possibility  during  the  entire 
history  of  the  system,  regardless  of  any  intermediate  states.  For  example,  in  the  "die  in  the 
funnel"  example,  we  may  be  sure  that,  whatever  happens,  the  die  will  not  come  out  the  top 
of  the  funnel.  FROB  is  able  to  predict  that  the  system  must  end  up  in  a  stable  state,  and 
therefore  it  must  traverse  a  path  in  the  state  transition  graph  which  ends  in  a  stable  state.  We 
believe  that  effective  qualitative  reasoning  requires  an  increased  use  and  development  of 
principles  of  this  sort,  and  a  decreased  reliance  on  simulation  techniques. 

A  natural  knowledge  engineering  approach  to  the  problem  would  be  to  write  down  rules 
that  simply  state  the  desired  prediction,  such  as  "A  small  object  released  inside  a  steep, 
large-mouthed  funnel  will  fall  out  the  bottom."  But  rules  of  this  kind  are  inadequate,  and 
have  rightly  been  rejected  by  previous  researchers.  Any  single  rule  of  this  kind  can  cover 
only  a  very  small  class  of  problems;  to  cover  large  classes  of  problems,  we  would  need  many 
separate  disconnected  rules.  In  particular,  a  rule  like  the  one  suggested  above  will  apply  only 
when  the  die  and  the  funnel  are  the  only  objects  involved.  As  soon  as  another  object  enters, 
the  conclusion  may  not  apply,  and  the  rule  gives  no  guidance.  In  technical  terms,  rules  of  this 
kind  are  not  compositional  across  objects  (Davis,  86b).  Even  when  dealing  with  only  a  die 
and  a  funnel,  if  we  allow  wide  ranges  in  the  shape  of  these  objects,  the  conclusion  will 
sometimes  apply  and  sometimes  not.  Since  there  is  no  simple,  general  rule  for  when  a  die 
will  come  out  the  bottom,  many  different  rules,  corresponding  to  special  geometric  cases, 
must  be  stated. 

There  are  three  problems  with  maintaining  a  knowledge  base  with  many  separate  rules 
corresponding  to  many  special  cases.  First,  in  the  long  run  it  will  lead  to  a  very  large 
knowledge  base,  and  thus  to  space  and  time  inefficiencies.  Second,  if  a  new  case,  not 
precisely  covered  in  pre-canned  categories,  comes  up,  the  system  cannot  even  begin  to  deal 
with  It.  Third,  this  approach  is  aesthetically  distasteful.  It  seems  clear  that  a  well-designed 
system  should  use  similarities  among  different  cases  of  a  die  falling  through  a  funnel,  and 
similarities  between  this  problem  and  similar  problems,  such  as  a  die  shot  through  a  tube,  or 
a  die  dropped  into  a  box.  These  analyses  of  these  cases  ought  to  have  more  in  common  than 
the  use  of  rules  which  are  syntactically  similar.  For  these  reasons,  we  have  rejected  the  use  of 
special  purpose  rules  of  this  kind  in  our  system. 

3.    Examples  and  Analysis 

We  propose  that  the  correct  analysis  of  the  "die  in  the  funnel"  example  takes  the 
following  form:  (i)  Due  to  the  topology  of  the  funnel,  if  the  die  goes  from  inside  it  to  outside 
it,  it  must  either  exit  the  top  or  exit  the  bottom,  (ii)  Since  the  die  is  dropped  from  rest  inside 
the  funnel,  it  cannot  have  the  energy  to  exit  the  top  of  the  funnel,  (iii)  There  is  no  stable 
resting  point  for  the  die  inside  the  funnel,  since  it  is  smaller  than  the  funnel's  mouth,  and  the 
funnel's  sides  are  steep,  (iv)  The  die  cannot  stay  forever  moving  within  the  funnel,  for  its 
kinetic  energy  will  eventually  be  dissipated.  Therefore,  the  die  must  exit  the  bottom  of  the 
funnel.  We  claim  that  in  most  cases  where  common  sense  predicts  that  the  die  will  come  out 
the  bottom  it  will  be  possible  to  carry  out  such  an  analysis,  and  to  support  the  substeps  by 
inferences  from  rules  of  fairly  broad  applicability.  (There  are  exceptions  where  this  analysis 
is  unnecessary  and/or  impossible;  for  instance,  when  the  die  is  dropped  from  a  great  height 
directly  above  the  mouth  of  the  funnel.)  Different  categories  of  problems  will  vary  in  the 
justifications  of  the  substeps. 

Related  problems  will  share  parts  of  the  analysis.  For  instance,  in  predicting  that  a  die 
in  a  small-necked  funnel  will  come  to  rest  at  the  top  of  the  neck,  we  may  use  the  identical 


arguments  (i)  that  the  die  must  either  exit  the  top.  exit  the  bottom,  or  stay  inside;  (ii)  thai  it 
cannot  exit  the  top;  and  (iv)  that  it  cannot  stav  inside  in  a  perpetual  state  of  motion.  The 
argument  (iii)  that  it  cannot  rest  stably  inside  the  funnel  must  be  modified  to  an  argument 
that  it  can  only  rest  stably  at  the  top  of  the  neck  of  the  funnel;  and  the  additional  argument 
must  be  made  that  it  cannot  exit  the  bottom  of  the  funnel,  since  the  orifice  is  too  small. 

This  analysis  avoids  both  problems  discussed  in  section  2.  We  can  avoid  analysing,  or 
even  determining,  the  states  of  motion  of  the  die  inside  the  tunnel;  all  we  need  to  determine 
is  that  the  die  cannot  rest  stably  inside.  Different  categories  of  problems  are  analysed  in 
similar,  but  not  identical  ways,  from  principles  of  broad  utility. 

In  the  remainder  of  this  section,  we  will  look  at  several  variations  of  the  "die  in  the 
funnel"  example,  and  show  how  an  analysis  of  this  sort  can  be  applied  to  each  of  them. 

We  begin  with  the  simple  case  shown  in  figure  4.  The  die  is  a  uniform  sphere.  The 
funnel  is  the  surface  of  revolution  about  a  vertical  axis  of  a  planar  figure  with  a  convex  inner 
side.  The  radius  of  the  die  is  less  than  the  radius  of  revolution  of  the  funnel.  The  steps  of  the 
argument  are  easily  filled  in.  (i)  The  top  and  bottom  of  the  funnel  are  the  only  orifices  of 
free  space  connecting  the  inside  of  the  funnel  with  its  outside.  Therefore,  if  the  die  is  to  go 
from  inside  to  outside,  it  must  go  through  the  top  or  the  bottom,  (ii)  Since  the  die  is 
spherical,  its  center  of  mass  is  in  its  interior.  Since  the  top  of  the  funnel  is  horizontal,  and 
directed  upward,  if  the  die  were  to  exit  it,  for  each  point  in  the  interior  of  the  die  there 
would  be  an  instant  of  time  when  it  was  above  the  top  of  the  funnel.  In  particular,  there 
would  be  an  instant  when  the  center  of  mass  would  be  above  the  top.  But  the  die  started  out 
from  rest  below  the  top  of  the  funnel,  and  there  is  no  source  of  additional  energy  for  the  die. 
Therefore,  the  die  cannot  come  out  the  top.  (iii)  By  a  simple  geometrical  argument,  the  die 
can  only  abut  the  inside  of  the  funnel  in  a  single  point.  A  rule  of  statics  asserts  that  a  uniform 
can  be  stably  supported  at  a  single  point  only  if  the  supporting  surface  is  horizontal  at  that 
point.  But  the  inner  surface  of  the  funnel  is  nowhere  horizontal.  Hence  there  is  no  resting 
place  for  the  die  inside  the  funnel.  The  full  details  of  this  argument  are  given  m  sections  6 
and  7. 

3.1.  Out  the  top,  out  the  bottom,  or  stay  inside 

We  now  consider  how  the  various  parts  of  this  argument  can  be  generalized  and 
modified.  Part  (i),  that  the  die  must  either  exit  the  top,  exit  the  bottom,  or  stay  inside,  does 
not,  of  course,  require  that  the  funnel  be  a  solid  of  revolution;  it  requires  only  that  the  funnel 
be  a  tube  with  only  two  orifices.  We  can  weaken  the  condition  still  further,  and  require  only 
that  all  orifices  other  than  the  top  or  the  bottom  be  too  small  to  let  the  die  go  through.  The 
determination  of  whether  a  die  can  go  through  a  hole  is  a  geometric  calculation  which  can  be 
easily  carried  out  for  a  variety  of  special  cases.  If  we  can,  by  other  means,  determine 
something  about  the  path  followed  by  the  die,  then  we  can  weaken  this  condition  still  further 
to  require  only  that  those  orifices  that  can  be  reached  by  the  die.  other  than  the  top  and 
bottom,  are  too  small.  For  example,  in  figure  5,  if  we  can  show  that  the  die  could  never 
reach  the  large  hole  in  the  side,  then  we  can  again  conclude  that  the  die  can  only  come  out 
the  top  and  the  bottom.  If  there  is  more  than  one  object  in  the  funnel,  then  the  rule  applies  to 
each  of  them  separately;  this  is  one  rule  which  is  compositional  across  objects  in  a  simple 
way. 

3.2.  Not  out  the  top 

Part  (ii),  the  argument  from  energy  conservation  that  the  die  cannot  come  out  the  top, 
depends  on  the  die  being  convex  and  on  the  center  of  mass  of  the  die  starting  out  below  any 
part  of  the  top.  Actually,  the  only  reason  we  need  convexity  is  to  establish  that  the  center  of 
mass  of  the  object  is  in  its  interior.  If  we  can  do  this  otherwise  —  for  example,  by  exact 
calculation,  or  by  establishing  that  the  shape  of  the  object  lies  within  some  sufficiently  small 
perturbation  of  a  convex  shape,  and  that  the  mass  is  not  all  concentrated  in  the  perturbation 
—  then  that  is  sufficient. 


We  can  weaken  the  requirement  on  the  shape  of  the  die  still  further  All  we  need  is 
that  the  ringed  filling  m  of  the  die  be  convex.'  The  ringed  filling  in  of  a  three-dimensional 
shape  5  is  defined  as  follows:  Consider  any  planar  cross  section  of  S  .  Let  C  be  any  simple 
closed  curve  that  lies  entirely  in  this  cross  section.  Let  p  he  a  point  in  the  plane  in  the  inside 
of  C .  Then  p  is  in  the  ringed  filling  in  of  S .  (Figure  6) 

Let  S  be  the  shape  of  some  object  O.  and  let  R  be  the  ringed  filling  m  of  S.  Assume  R 
is  convex.  Clearly,  R  is  equal  to  the  convex  hull  of  S ,  so  R  contains  the  center  of  mass  of  0. 
Let  C  be  a  closed  curve  lying  in  S  and  in  a  plane  containing  the  center  of  mass  of  O.  It  can 
easily  be  seen  that,  if  C  goes  through  a  planar  surface,  then  so  does  every  point  inside  C. 
Thus,  if  O  exits  the  top  of  the  funnel,  then  the  center  of  mass  of  O  must  likewise,  and  the 
proof  goes  through.  This  enables  us  to  establish  step  (ii)  for  such  shapes  as  a  torus,  or  a 
wiffle  ball,  or  a  cratered  convex  shape. 

We  can  also  apply  the  same  argument  in  cases  where  we  have  two  or  more  dice  bound 
together.  If  the  union  of  the  shapes  of  the  dice  is  always  convex,  or  its  ringed  filling  in  is 
always  convex,  then  the  same  argument  applied  to  the  center  of  mass  of  the  pair  will  go 
through  in  very  nearly  the  same  way.  For  instance,  if  one  of  the  dice  contains  the  other,  and 
has  a  convex  exterior,  as  in  figure  3,  then  it  will  not  come  out  the  top. 

However,  the  argument  cannot  be  extended  in  any  straightforward  way  to  handle 
boomerang  shaped  objects  (Figure  7).  The  problem  is  that  the  center  of  mass  of  a  boomerang 
is  inside  its  curve,  and  that  a  boomerang  actually  can  come  out  of  the  top  of  funnel  while 
keeping  the  center  of  mass  always  lower  than  the  top.  (This  is  the  principle  used  by  high 
jumpers.)  Nor  does  it  apply  if  there  are  several  dice  which  are  not  bound  together.  In  figure 
8,  for  example,  it  is  possible  that  die  B  will  eventually  collide  with  die  A  in  such  a  way  that 
A  gets  enough  energy  to  escape  out  the  top.  How  can  we  handle  these  cases? 

One  approach  might  simply  be  to  exclude  these  behaviors  by  fiat  on  the  grounds  that 
they  would  be  very  unlikely.  I  am  not  entirely  happy  with  this  approach,  not  just  because  the 
rules  would  not  be  strictly  true,  but,  more  importantly,  because  it's  not  clear  to  me  that  these 
behaviors  always  are  unlikely.  In  particular,  if  you  did  an  experiment  dropping  elastic  dice.t 
such  as  Ping  Pong  balls,  inside  a  large  funnel,  my  guess  would  be  that  one  of  them  would 
come  out  the  top  fairly  often.  .My  instincts  about  elastic  boomerangs  are  much  feebler,  but  I 
wouldn't  bet  much  money  against  it  ever  coming  out  the  top.  In  short,  I  think  that  the 
intuitive  prediction  in  figures  7  and  8  that  A  comes  out  the  bottom  rests  on  an  assumption 
that  the  dice,  like  most  objects,  are  fairly  inelastic. 

This  being  the  case,  we  should  postulate  that  the  dice  have  a  coefficient  of  elasticity  of 
^.<1  and  build  an  argument  around  that.  This  turns  out  to  be  trickier  than  one  might 
imagine;  inelastic  collisions  are  rather  strange.  On  first  consideration,  you  might  think  that 
an  argument  can  be  constructed  along  the  following  lines:  Let  us  postulate  that  the 
boomerang  is  initially  dropped  over  a  section  of  the  funnel  which  is  substantially  below  the 
top  of  the  funnel.  To  be  precise,  let  d  be  the  diameter  of  the  boomerang,  and  let  h  be  the 
distance  from  the  top  of  the  funnel  to  the  point  of  the  funnel  underneath  the  boomerang. 

Then  assume  h> .  Now  the  following  argument  seems  plausible:  When  the  boomerang 

strikes  the  side  of  the  funnel,  it  will  have  at  most  energy  E  =  mgh.  (m  =  mass  of  boomerang; 
g  =  acceleration  of  gravity.)  After  the  collision  with  the  funnel,  it  will  have  at  most  energy 
E'  =  \Lmgh.  At  the  moment  of  collision,  the  height  of  the  center  of  mass  of  the  boomerang 
above  the  funnel  is  at  most  d.  Hence,  the  maximum  height  attainable  bv  the  center  of  mass 

E' 

after  the  collision  is  at  most  d+ =  ^+  ^./i.  But  d-^  \x.h<h- d,  so  that  after  the  collision,  the 

mg 

boomerang  cannot  get  high  enough  for  any  part  of  it  to  reach  the  top. 


*  I  made  up  this  terminology.  Please  let  me  know  if  there  exisis  a  more  standard  terra  for  this. 
"  By    elastic'  I  mean  that  they  lose  little  energy  in  collisions,  not  that  they  are  flexible 


The  argument  could  be  made  tighter,  and  the  condition  on  h  made  considerably  weaker, 
but  that  IS  a  minor  point.  In  fact,  the  argument  is  intrinsically  flawed.  It  is  not  the  case  that  a 
collision  of  an  inelastic  object  with  a  fixed  object  necessarily  results  in  a  loss  of  1  -  p.  of  its 
incoming  kinetic  energy.  It  may,  in  fact,  result  in  no  loss  of  energy  whatever.  Consider  the 
situation  shown  in  figure  9.  No  impact  takes  place;  hence,  no  energy  is  lost  during  impact. 
No  energy  is  lost  to  the  normal  forces,  since  normal  forces  are  non-dissipative.  Hence,  in  the 
absence  of  friction,  no  energy  whatever  is  lost,  and  the  ball  will  fly  from  the  wall  at  the  same 
velocity  at  which  it  struck  the  wall.  If  the  ball  comes  in  with  the  correct  forespin,  no  energy 
will  be  lost  in  the  presence  of  friction,  either.  Other  anomalous  cases  which  must  be 
considered  in  formulating  a  general,  useful  rule  of  collisions  over  extended  time  include 
figure  10,  in  which  one  inelastic  object  sideswipes  another,  with  very  little  energy  loss;  and 
figure  11,  in  which  a  ball  thrown  up  a  slanting  rooftop  comes  down  again  at  very  nearly  the 
same  speed. 

I  think  that  it  probably  possible  to  formulate  true  rules  which  put  lower  bounds  on  the 
energy  loss  during  a  collision  over  extended  time.  I  would  guess,  for  example,  that  a  rule 
something  like  the  following  is  true:  Let  object  A  be  moving  with  velocity  v  towards  fixed 
object  B.  Let  the  normal  to  the  surface  of  S  at  every  point  where  A  comes  into  contact  of  B 
be  within  some  small  angle  <i)  of  -v.  Let  the  coefficient  of  elasticity  of  A  be  some  small 
value  p..  Let  t.  be  the  moment  when  A  comes  into  contact  with  B,  and  let  r,  be  any  moment 
after  A  has  attained  a  translational  velocity  in  a  direction  opposed  to  v  Let  E  be  the  kinetic 
energy  associated  with  the  translational  motion  of  A  at  f, ;  and  let  £"  be  the  kinetic  energy  lost 
to  conservative  forces  (gravity)  between  r,  and  r,.  Assume  that  A  does  not  come  into  contact 
with  any  object  other  than  B  between  f,  and  r,.  Then  the  collision  will  dissipate  some  fraction 
of  the  energy  E-E' .  This  fraction  is  bounded  below  by  some  decreasing  function  of  p.  and 
<().  (Figure  12).  In  other  words,  head-on  collisions  lead  to  a  loss  of  the  energy  of  translational 
motion  (rotational  energy  may  survive  intact.)  As  I  say,  this  seems  plausible  and  is  probably 
true. 

For  the  purposes  of  forming  a  logical  system,  it  is  not  necessary  to  specify  the  exact 
function  which  relates  \i.  and  4)  to  the  energy  loss  above.  We  can  just  give  the  function  a 
name,  and  use  it  in  our  problem  statements  and  proofs.  In  an  actual  system,  of  course,  it 
would  be  important  to  say  something  quantitative  about  the  function.  It  is  probably  not 
terribly  important  how  the  system  designer  finds  this  quantitative  information.  He  could 
work  it  out  analytically,  or  derive  it  experimentally,  or  guess  at  it,  or  have  the  program 
induce  it  from  observations.  It  should  not  be  hard  to  come  up  with  some  values  that  are 
conservative  enough  to  be  safe,  and  liberal  enough  to  be  useful. 

Unfortunately,  this  rule  is  no  use  in  the  funnel  example,  since  the  sides  of  the  funnel 
may  well  be  nearly  parallel  to  the  direction  of  incoming  motion  (downwards).  In  fact,  I  have 
not  been  able  to  find  any  adequate  analysis  of  this  example  or  of  the  example  with  rwo  dice 
dropped  inside  a  funnel,  or  of  the  very  basic  problem  of  bounding  the  maximum  horizontal 
translation  possible  when  an  eraser  is  dropped  on  the  floor.  The  analysis  of  collisions  over 
extended  time  is  the  most  important  gap  in  our  logical  theory. 

3.3.    No  resting  point  inside 

Part  (iii)  of  our  argument,  that  the  die  cannot  rest  inside  the  funnel,  depended  in  our 
first  example  on  the  very  strong  assertions  that  the  die  was  a  uniform  sphere  and  that  it  could 
contact  the  funnel  only  in  a  single  point.  We  can  easily  achieve  a  slight  generalization  to 
nearly  uniform,  nearly  spherical  dice,  which  behave  similarly,  not  surprisingly.  One  can 
establish  the  following  formula:  Let  6  be  the  slope  of  the  support;  let  p.  be  the  coefficient  of 
friction;  let  i>  be  the  maximum  angle  berween  the  line  from  the  center  of  mass  to  a  point  on 
the  surface  and  the  normal  to  the  surface  at  that  point.  (Figure  14)  The  ball  can  stand  still 
only  if  p.atan(9)  and  (l)>e.  We  may  further  easily  generalize  the  result  by  replacing  the 
condition  that  the  objects  contact  in  a  single  point  with  the  weaker  condition  that  the  surface 
normal  at  all  contact  points  are  parallel.  Thus,  we  can  use  these  formulas  to  predict  that  a 


regular  icosahedral  die  will  roll  or  slide  down  a  funnel. 

We  can  similarly  handle  the  case  of  a  die  which  is  a  spherical  shell  and  contains  another 
die  inside.  (Figure  15)  Here  things  are  made  easier  by  the  fact  that  inner  die  will  tend  to  slide 
to  the  bottom  of  the  outer  die,  and  exert  a  torque  toward  rolling.  The  system  can  rest 
staticailv  only  if  there  is  a  position  in  which  joint  center  of  mass  of  the  two  dice  is  located 
directly  above  the  contact  point  of  the  outer  die  with  its  support,  and  in  which  the  inner  die 
can  rest  stably  inside  the  outer  die.  It  is  not  clear  whether  the  conditions  which  determine  this 
should  be  thought  of  as  high  level  axioms,  or  whether  they  should  be  derived  from  a  force 
analysis.  Either  is  possible. 

Once  we  get  to  dice  which  can  contact  the  funnel  in  several  points  with  different  surface 
normals,  the  analysis  becomes  sufficiently  harder  that  it  seems  reasonable  to  forget  about 
rolling,  and  just  look  for  conditions  sufficient  to  produce  sliding.  The  wider  the  range  of  the 
horizontal  component  of  the  surface  normals  at  contact  points,  the  steeper  the  slope  must  be, 
for  the  normal  forces  at  the  various  contact  points  will  tend  to  act  against  each  other;  hence, 
they  will  be  larger  than  they  would  be  otherwise;  hence,  the  friction  forces  will  be  larger. 
The  following  rule  holds:  Let  A  be  in  contact  with  B.  Let  9  be  the  minimum  slope  of  the 
surface  of  fl  at  a  contact  point.  Consider  the  horizontal  components  of  the  surface  normals  of 
B  at  the  contact  points,  and  assume  that  there  is  some  direction  which  lies  within  some  small 
angle  <i>  of  all  these  horizontal  components.  (Figure  16)  Let  the  coefficient  of  friction  be  p.. 
If  jjL<cos<t)tand,  then  A  will  slide  down  S.  This  formula  can  be  rigorously  shown  by  a  simple 
extension  of  the  usual  proof  that  a  block  on  an  inclined  plane  slides  if  ^i,<tane. 

The  geometric  condition  that  the  surface  normals  at  the  contact  points  are  all  nearly 
parallel  in  their  horizontal  components  can  be  established  in  various  ways.  For  example,  one 
can  get  this  by  bounding  the  diameter  of  the  die,  and  the  curvature  of  inner  surface  of  the 
funnel;  or  by  special  formula  for  some  class  of  shape  description;  etc. 

The  above  formula  will  work  just  as  well  for  multiple  dice,  if  we  consider  all  the  contact 
points  with  the  funnel  of  some  collection  of  dice  connected  by  abutments.  Thus,  we  can  easily 
show  that  if  a  collection  of  dice  is  so  small  that,  strung  end  to  end,  they  can  only  reach  points 
on  the  funnel  with  surface  normals  with  close  horizontal  components,  then  they  must  fall 
through.  It  IS  undoubtedly  possible  to  establish  stronger  rules  with  a  little  more  work. 

By  combining  all  the  different  ways  in  which  the  results  (i),  (ii),  and  (iii)  may  be 
established  together,  and  all  the  ways  in  which  their  geometrical  preconditions  may  be 
satisfied,  we  get  a  rich,  interconnected  body  of  results,  all  with  the  conclusion,  "The  die  falls 
out  the  bottom  of  the  funnel." 

4.    The  Underlying  Logic 

We  can  now  address  the  main  question  of  this  paper:  What  kind  of  logical  framework  is 
necessary  to  carry  out  the  reasoning  described  in  the  previous  section?  We  divide  our 
discussion  into  geometry,  temporal  logic,  and  physics.  We  end  this  section  by  discussing  the 
fact  that  Newtonian  physics  is  not  fully  deterministic,  a  fact  with  some  implications  for  the 
logic  of  the  domain. 

4.1.    Geometry 

The  geometry  of  the  world  is  taken  to  be  R^.  The  space  occupied  by  an  object  in  some 
particular  position  is  a  subset  of  /?';  it  is  assumed  to  be  closed,  connected,  and  normal.  By 
"normal",  we  mean  that  no  part  of  the  object  is  a  two  dimensional  lamina,  or  a  one 
dimensional  curve;  technically,  the  shape  of  the  object  is  equal  to  the  closure  of  its  interior. 
The  fact  that  objects  are  closed  means  that  abutting  objects  actually  share  some  boundary 
points;  though  this  seems  a  little  strange,  it  does  not  cause  any  actual  problems  in  the  logic. 
(See  (Fleck,  in  preparation)  for  an  interesting  topology  that  avoids  this  anomaly.)  We  will 
further  assume  that  the  boundaries  of  objects  are  everywhere  smooth,  with  a  well-defined 
tangent  plane.   It  is  probably  possible   to  do  without  this  assumption,  but  it  simplifies  the 


logic. 

The  arguments  in  section  3  used  several  different  kinds  of  geometric  reasoning.  These 
included  the  following: 

•  The  ability  to  name  and  describe  particular  point  sets  that  are  connected  to  objects  and 
are  useful  for  physical  reasoning,  such  as  the  top  and  bottom  of  a  funnel.  These  are 
called  "pseudo-objects"  in  our  system;  the  problem  of  constructing  them  effectively  is 
the  same  as  the  problem  of  constructing  metric  diagrams  in  FROB  (Forbus,  79). 

•  The  use  of  topological  predicates.  For  example,  the  funnel  forms  topologically  a  box 
with  two  orifices,  and  the  die  starts  out  inside  the  box.  (Predicates  of  this  kind  are 
examined  in  detail  in  (Davis,  84b).) 

•  The  use  of  a  property  quantified  over  all  irregularities  of  a  certain  kind  in  an  object. 
For  instance,  we  may  wish  to  say  that  the  funnel  has  no  holes  large  enough  to  let  the  die 
through  other  than  the  top  and  bottom. 

•  The  use  of  special  shapes,  such  as  spheres  and  surfaces  of  rotation,  and  the  use  of  rules 
peculiar  to  these  shapes. 

•  The  use  of  metric  dimensions,  and  the  statement  of  inequalities  on  them.  For  example, 
we  wish  to  say  that  the  radius  of  the  die  is  less  than  the  radius  of  the  funnel. 

•  The  bounding  of  the  range  of  the  surface  normal  over  a  part  of  an  object's  surface.  For 
example,  we  wish  to  say  that  the  slope  of  the  funnel  is  everywhere  positive  in  its  inner 
surface;  or  that,  for  each  region  of  a  specified  size  in  the  inner  surface,  there  is  a 
central  direction  such  that  the  horizontal  component  of  the  surface  normal  lies  within 
some  small  angle  of  the  central  direction 

•  The  use  of  convexity  and  related  properties. 

Thus,  any  adequate  geometric  language  will  have  to  be  strong  enough  that  these,  or 
most  of  these,  can  be  either  expressed  directly  or  inferred.  One  could,  of  course,  simply  use 
a  language  of  shape  which  allowed  all  these  kinds  of  information  to  be  stated  directly  and  . 
conjoined.  Such  an  approach  does  not  seem  to  be  on  the  right  track;  it  would  be  very  hard,  in 
such  a  language,  to  do  things  like  checking  a  shape  description  for  consistency,  or  comparing 
two  descriptions  for  possible  identity.  I  would  rather  find  a  more  systematic,  constructive, 
and  restricted  language,  and  I  am  looking  into  languages  which  combine  some  kind  of 
constructive  solid  geometry  with  facilities  for  describing  surface  characteristics.  The  real 
problem  is  to  make  the  language  flexible  enough  to  accomodate  the  desired  kinds  of 
uncertain  and  qualitative  description.  For  instance,  it  is  clear  that  the  concepts  of  "grain-size" 
or  tolerance"  used  in  (Davis,  86a)  and  in  (Requicha,  83)  are  not  adequate  in  this  domain, 
since  they  leave  the  surface  orientation  completely  unconstrained. 

4.2.    Temporal  Logic 

Our  temporal  logic  follows  that  of  (McDermott,  82)  fairly  closely.  A  scene  is  an 
instantaneous  snapshot  of  the  universe.  In  our  restricted  domain,  a  scene  consists  precisely  of 
a  specification  of  the  positions  and  velocities  of  all  objects.  A  chronicle  is  function  from  the 
time  line  to  scenes."  The  time  line  is  taken  to  be  R^. 

Unlike  (McDermott,  82),  we  do  not  need  the  concept  of  branching  time.  However,  it 
may  be  useful  to  use  multiple  chronicles  for  expressing  physical  principles  in  terms  of 
counterfactuals  (e.g.  "A  bullet  in  free  fall  is  always  lower  than  it  would  be  if  it  went  in  pure 
linear  motion ')  Note  that  the  use  of  the  counterfactual  means  that  chronicles  must  range  over 


"  Termmolcgical  note:  What  I  call  a  "scene"  here  and  in  my  previous  work  (Davis.  84a)  was  called  a 
"state'  in  (McDermott.  82),  though  "state"  was  used  with  a  different  meaning  in  (Chainiak  and 
McDermon.  85).  It  is  prettv-  much  the  same  thing  as  the  "situation'  of  (McCarthy  and  Hayes,  69).  but 
embedded  in  a  continuous,  rather  than  a  discrete,  time-line.  My  use  of  "chronicle"  accords  with 
(.McDermon.  82).  In  mv  own  previous  work  (Davis.  84a),  I  called  this  a  "history  ".  but  I  have  repented 
of  this  usage,  since  Pat  Hayes  (1978a)  long  ago  used  "history  '  to  means  something  quite  dilferent. 
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physically  impossible  chains  of  events,  as  well  as  physically  possible  ones.  (See  (Davis,  84a) 
and  (Shoham,  86)  for  further  examples.) 

In  the  particular  examples  of  section  3,  the  only  temporal  predicates  needed  relate 
purely  to  order  relationships.  For  other  predictions,  such  as  determining  roughly  when  the 
die  may  be  expected  to  come  out  the  bottom  of  the  funnel  (in  a  second?  in  a  minute?  next 
year?)  actual  measurements  on  time  must  be  used. 

It  may  be  noted  that  the  "frame"  or  "persistence"  problem  of  determining  what  remains 
true  over  time  ((McCarthy  and  Hayes,  69),  (McDermott,  82),  (Shoham,  86))  requires  no 
special  treatment  in  our  logic.  We  can  avoid  this  problem,  not  by  virtue  of  any  special 
cleverness  on  our  part,  but  by  virtue  of  the  structure  of  the  domain.  The  predicates  in  the 
domain  are  divided  into  two  classes.  The  first  class  includes  predicates  which  depend  on 
position  and  velocity  of  objects.  These  are  not  assumed  to  remain  constant  over  any  interval 
unless  they  can  be  proven  to  be  so.  The  second  class  includes  structural  predicates,  depending 
only  on  the  shape  and  other  material  properties  of  the  objects.  These  are  assumed  to  be 
constant  over  the  problem,  and  so  are  defined  atemporally.  Similar  considerations  would 
seem  to  apply  to  any  closed  world,  complete  physical  theory;  it  is  not  clear  why  the  frame 
problem  should  ever  create  trouble  in  such  a  domain.  (The  closed  world  assumption  is 
explicitly  introduced  in  our  theory  through  the  predicate  "isolated(00,C)",  which  asserts  that, 
during  C,  no  mobile  object  in  the  set  of  objects  00  ever  comes  into  contact  with  any  object 
outside  00.  See  section  6.3,  axiom  6.3.1.6.) 

4.3.    Physics 

Our  worlds  consist  of  a  finite  set  of  solid  objects  which  move  about  in  space  through 
time.  Objects  are  completely  rigid,  and  indestructible;  the  set  of  existent  objects  always  stays 
the  same.  The  interior  of  objects  may  not  overlap. 

Objects  are  subject  to  four  kinds  of  forces;  a  uniform  downward  gravitational  force; 
normal  forces,  which  act  to  prevent  objects  from  overlapping;  friction;  and  a  weak  drag 
force,  which  serves  to  dissipate  energy  of  motion  when  no  other  force  is  available.  The  drag 
force  is  assumed  to  be  very  much  smaller  than  any  other  force.  We  need  to  assume  its 
existence  in  order  to  postulate  that  no  system  can  persist  indefinitely  in  a  state  of  motion. 
Friction  by  itself  will  not  do  this  for  us;  a  rolling  ball  loses  no  energy  to  friction. 

Certain  objects  are  fixed.  A  fixed  object  does  not  move,  no  matter  what  forces  are 
applied  to  it.  One  particular  fixed  object,  which  is  part  of  every  chronicle,  is  the  ground.  The 
ground  is  a  fixed  object  of  infinite  horizontal  extent;  our  model  is  weak  on  astronomy, 
though  strong  on  physics.  Any  vertical  line  intersects  the  surface  of  the  ground  in  a  unique 
point;  all  points  below  the  surface  are  part  of  the  ground. 

Objects  have  two  other  internal  properties  besides  their  shape:  a  distribution  of  mass 
within  their  interior,  and  a  coefficient  of  elasticity,  which  determines  how  the  object  behaves 
in  a  collision.  Any  pair  of  objects  have  a  coefficient  of  friction,  which  determines  the  frictive 
forces  when  the  objects  are  in  contact. 

Thus  a  particular  world  is  determined  by  enumerating  the  objects  in  the  world  with  their 
shapes,  mass  distributions,  coefficients  of  restitution,  and  pairwise  coefficients  of  friction. 
When  a  collision  takes  place,  the  velocity  is  discontinuous.  By  convention,  we  will  define  the 
velocity  of  an  object  in  a  collision  scene  as  the  limit  of  its  velocity  from  previous  times.  (I 
assume,  without  proof,  that  this  limit  always  exists.)  Thus,  we  may  specify  as  the  initial 
conditions  of  a  problem  that  a  collision  is  taking  place. 

We  need  a  fairly  flexible  ontology  for  forces.  The  problem  is  that  one  object  may  exert 
force  on  another  at  a  single  point,  along  a  curve,  or  over  a  surface  (figure  17).  In  simple 
physics  problem,  one  generally  considers  the  resolvant  force  as  acting  at  a  single  point,  but 
this  is  not  feasible  in  the  more  general  case,  since  the  resolvant  force  may  not  be 
perpendicular  to  the  surface  there.  The  solution  is  to  consider  a  force  on  an  object  as  a 
distribution  of  force  over  the  volume  of  the  object.    That  is,  given  any  subset  of  the  object 
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shape,  a  certain  force  is  applied  to  the  object  over  that  subset.  (If  we  wished  to  treat 
collisions  in  terms  of  forces  as  well,  then  we  would  need  to  treat  forces  as  distributions  of 
momentum  transfer  over  space  and  time.  However,  my  guess  is  that  it  is  better  to  treat  the 
momentum  transfer  in  collisions  separately  from  ordinary  forces.) 

The  physical  deductions  which  must  be  made  in  the  examples  of  section  3  include  the 
following: 

•  Determining  whether  a  set  of  objects  can  attain  a  stable  scene  while  certain  geometric 
conditions  hold. 

Finding  constraints  on  the  location  of  the  center  of  mass  of  an  object  or  a  set  of  objects. 
Resolving  a  set  of  forces,  and  determining  motion  under  those  forces. 
Predicting  that  the  existing  structure  of  contacts  between  objects  will  change. 
Determining  that  a  collision  will  take  place. 
Predicting  the  effect  of  a  collision. 

Determining  whether  a  chronicle  violates  a  conservation  law,  particularly  conservation 
of  energy. 

•  Characterizing  the  paths  that  an  object  can  take  without  coming  to  overlap  other  objects. 

4.4.    Non-determinism 

Contrary  to  popular  impression,*  Newtonian  physics  is  not  deterministic.  It  is  not  the 
case  that,  given  a  starting  state,  all  subsequent  history  is  determined;  it  is  nearly  the  case,  but 
there  are  some  frustrating  exceptions.  These  exceptions  arise  from  two  facts:  first,  the 
distribution  of  forces  in  a  given  situation  may  be  ill-defined;  second,  what  happens  at  a 
collision  may  be  ill-defined.  We  will  discuss  each  of  these  briefly. 

If  a  stationary  object  is  supported  at  more  than  three  points  (more  than  two,  in  two 
dimensions),  then  the  distribution  of  forces  among  these  support  points  may  be  ill  defined.  In 
the  two  dimensional  case,  there  are  only  two  constraints  —  that  the  support  forces  balance 
the  weight  of  the  object,  and  that  their  torques  sum  to  zero  —  and  these  are  not  sufficient  to 
fix  the  magnitudes  of  the  three  forces  involved.  For  example,  in  figure  18,  the  block  X  may 
be  puning  all  its  weight  on  B;  or  half  its  weight  on  A  and  half  on  C;  or  one-third  on  each;  or 
any  other  symmetric  distribution.  There  is  no  way  to  tell. 

Ordinarily,  it  makes  no  difference;  however  the  forces  are  distributed,  block  X  will  stay 
put.  Under  some  special  circumstances,  however,  it  may  make  a  different,  because  it  affects 
the  friction  force.  Imagine,  for  example,  that  the  situation  of  figure  18  is  tilted  toward  the 
line  of  sight.  Assume  that  blocks  A  and  C  are  fixed  to  the  table  and  that  the  coefficient  of 
friction  between  blocks  A,  B,  C.  and  X  is  high,  but  the  coefficient  of  friction  between  the 
blocks  and  the  table  is  low.  Now  the  distribution  of  forces  will  make  a  difference.  If  X  has  all 
of  its  weight  on  A  and  C,  it  will  stay  put,  and  B  will  slide  out  from  under.  If  X  has  all  its 
weight  on  B,  then  X  and  B  will  slide  together.  Intermediate  distributions  may  have  other 
effects;  for  example,  for  some  distributions  and  coefficients,  A  and  C  will  hold  X  in  place, 
and  X  will  hold  B. 

Similar  things  can  happen  if  one  object  is  "pinched"  by  another,  or  by  two  others  which 
are  fixed  with  respect  to  one  another.  In  figure  19,  for  example,  there  is  no  way  to  know 
what  are  the  normal  forces  between  A  and  B  on  the  sides;  all  that  can  be  known  is  that  they 
balance.  Again,  it  may  make  a  difference  in  behavior  through  the  frictive  forces;  if  A  is 
jammed  into  B,  and  the  normal  forces  are  large,  then  the  friction  may  be  sufficient  to  hold  A 
up.  Newtonian  physics  does  not  predict  how  tightly  a  peg  is  held  inside  a  hole.  Finally,  if  a 
coffee  can  is  spun  around  its  vertical  axis  on  a  table  top,  then  the  radial  distribution  of  the 
normal  force  will  affect  the  torque  exerted  by  the  friction,  and  consequently  will  affect  the 


At  least.  It  used  to  be  my  impression. 


time  needed  for  the  can  to  stop.    If  the  can  is  resting  on  its  center,  it  will  spin  a  long  time;  if 
It  IS  resting  on  its  edges,  it  will  stop  soon. 

The  other  source  of  ambiguity  in  Newtonian  physics  occurs  in  collisions.  A  collision 
involving  more  than  two  objects  is  almost  always  underdetermined.  The  only  constraints  are 
conservation  of  energy  (relative  to  the  coefficients  of  elasticity),  conservation  of  momentum 
and  angular  momentum,  and  the  constraint  that  the  direction  of  momentum  transfer  be 
normal  to  the  surface  of  collision  (ignoring  friction).  In  general,  these  are  not  sufficient  to 
determine  the  velocities  following  the  collision.  (See  figure  20.)  It  may  be  observed  that  a 
simultaneous  collision  among  several  objects  may  be  viewed  in  several  ways  as  the  limiting 
case  of  verv  closely  timed  two  object  collisions,  and  that  the  behavior  in  these  various  cases 
do  not  have  the  same  limit.  For  example,  in  figure  21,  the  behavior  of  the  system  if  block  A 
strikes  B  immediately  before  C  strikes  B  is  quite  different  from  the  behavior  if  C  strike  B 
before  A.  The  well-known  toy  with  a  row  of  steel  balls  swinging  as  pendulums  from  strings 
is  quite  misleading,  in  this  respect;  it  happens  that  sensitivity  to  timing  goes  away  in  the 
linear  motion  of  objects  of  equal  mass. 

A  final  underdetermined  case  is  where  two  objects  collide  along  a  non-planar  surface. 
In  this  case,  the  ambiguity  in  force  distribution  results  in  an  ambiguity  in  the  direction  of 
momentum  transfer,  and  hence  in  an  ambiguity  in  out-going  velocities.  (Figure  22) 

If  Newtonian  physics  does  not  guarantee  a  unique  solution,  can  we  be  sure  that  it 
guarantees  the  existence  of  a  solution?  I  do  not  know  whether  this  problem  has  been  solved 
in  general.  It  is  certainly  not  the  case  that  any  starting  scene  can  be  extrapolated  to  a 
chronicle.  Consider,  for  example,  a  scene  where  object  A  is  tightly  enveloped  in  object  B, 
where  the  two  objects  have  initially  different  velocities,  and  both  are  perfectly  elastic  (figure 
23);  Then  there  is  no  solution;  there  is  no  behavior  which  conserves  both  energy  (required 
by  the  perfect  elasticity)  and  momentum.  Of  course,  there  is  something  funny  about  this 
example;  namely,  there  is  no  way  to  get  into  this  state.  I  have  not  been  able  to  construct  an 
example  that  forces  itself  into  such  a  situation,  nor  to  prove  that  there  does  not  exist  any  such 
example.  I  suspect  that  no  such  example  exists,  and  that  the  Newtonian  theory  of  smooth 
solid  objects  never  gets  into  trouble.  If  there  were  such  an  example,  it  would  be  interesting;  it 
would  be  a  Gedanken  experiment  proving  that  objects  cannot  always  be  perfectly  rigid  and 
elastic. 

Since  the  ambiguity  of  Newtonian  physics  cannot  be  avoided,  we  have  to  accept  it.  It 
has  the  effect  of  complicating  some  predicates.  For  example,  we  cannot  have  a  single 
predicate  which  asserts  that  a  scene  is  stable.  Instead  we  need  two;  one  which  asserts  that  a 
scene  is  possibly  stable,  and  one  which  asserts  that  it  is  necessarily  stable.  A  block  sitting 
motionless  on  a  horizontal  table  is  necessarily  stable;  the  peg  in  figure  19  is  possibly  stable. 

A  hard  question  arises  as  regards  the  scene  in  figure  19.  We  have  observed  that  it  is 
indeterminate  whether  the  peg  stays  in  or  falls  out.  However,  if  the  peg  has  already  been 
motionless  for  a  while,  is  it  be  reasonable  to  posit  that  it  will  remain  motionless?  or  is  it 
possible  that  the  peg  may  stay  still  a  while,  then  slide  down  some,  then  stay  still  some  more, 
and  so  on''  If  we  posit  that  it  remains  motionless,  then  we  have  a  circumstance  where,  though 
the  scene  does  not  uniquely  determine  the  subsequent  chronicle,  chronicles  over  small 
intervals  may  do  so.  The  answer  depends,  at  least  in  part,  on  what  view  we  take  of  physical 
properties  that  lie  outside  the  range  of  our  theory,  such  as  temperature;  for  the  true  behavior 
of  the  peg  is  governed  by  these  properties.  If  these  properties  are  required  to  stay  constant 
throughout  a  chronicle,  then  it  is  reasonable  to  extrapolate  future  stability  from  present 
stability.  If  they  are  allowed  to  vary,  then  it  is  not  reasonable.  For  example,  the  peg  may 
move  down  its  hole  in  starts  and  stops  if  the  temperature  fluctuates  rapidly;  it  will  move  in 
low  temperature,  and  stop  in  high  temperature.  On  the  whole,  I  think  it  more  reasonable  to 
allow  this  kind  of  behavior  than  to  exclude  it. 
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5.    Ontology 

The  remainder  of  this  paper  explicitly  constructs  a  first-order  language  and  a  set  of 
axioms  adequate  to  analyse  the  first  and  simplest  example  of  section  3.  Section  5  constructs 
the  formal  ontology  needed  for  solid  object  physics.  Here,  we  go  somewhat  beyond  the 
minimal  requirements  of  our  chosen  example,  and  discuss  concepts  such  as  force  which  are 
necessary  only  in  more  complex  examples.  Section  6  defines  the  non-logical  terms  of  the 
language  and  states  the  axioms  needed  to  analyse  the  particular  example.  Section  7  carries 
out  this  analysis.  I  have  not  carried  out  an  analysis  at  this  level  of  detail  for  any  of  the  other 
examples  in  section  3.  I  believe  it  can  be  done,  but  it  is  easy  to  be  over-optimistic  on  such 
points. 

Only  the  reader  with  the  highest  degree  of  technical  interest  will  wish  to  work  through 
these  three  sections  in  detail.  The  more  casual  reader,  however,  might  like  to  get  a  feel  for 
what  I  am  doing  by  skimming  a  few  of  the  definitions  and  axioms  in  section  5  and  each 
subsection  of  section  6,  and,  particularly,  by  reading  the  exact  problem  statement  at  the 
beginning  of  section  7. 

Our  physical  ontology  contains  numerous  types  of  elements:  instants  of  time,  quantities 
of  mass,  quantities  of  energy,  points,  vectors,  vector  fields,  mappings,  movements,  objects, 
scenes,  pseudo-objects,  forces,  chronicles,  and  sets  of  these.  Unlike  section  3,  we  will  assume 
in  these  sections  that  all  objects  have  uniform  mass  distributions.  It  would  be  easy  to  extend 
the  logic  so  this  assumption  can  be  made  or  dropped  explicitly. 

Time,  mass,  and  energy  are  all  measure  spaces,  which  we  may  view  as  real  numbers.  It 
is  not  necessary  to  consider  branching  time  in  our  domain.  The  value  "=<:"  (infinity)  is  an 
acceptable  value  for  time. 

We  use  R^  as  our  underlying  geometry.  Points  and  vectors  have  their  usual  meaning. 
We  will  use  the  phrase  "geometric  entity"  to  mean  a  point,  a  set  of  points,  a  vector,  or  a 
vector  field.  We  specify  changes  in  positions  using  rigid  mappings;  these  are  adequate,  since 
we  are  dealing  only  with  rigid  objects.  A  rigid  mapping  is  the  composition  of  a  translation 
and  a  rotation;  reflections  are  not  allowed.  A  general  velocity  is  a  derivative  of  a  rigid 
mapping.  A  general  velocity  is  the  composition  of  a  linear  velocity  and  an  angular  velocity 
around  a  specified  axis. 

Definition  5.1:  A  point  set  is  a  subset  of  R^. 

Definition  5.2:  A  vector  field  is  a  pair  <XX.f>,  where  XX  is  a  point  set  and  /  is  a 
function  from  XX  onto  vectors.  We  use  the  function  "base(P)"  and  "■valueiP ,X)"  to  recover 
the  components  of  the  vector  field.  If  P  =<XXJ>  is  a  vector  field,  then  base(F)  =  XX, 
and,  for  XiXX,  value(F,X)  =  f(X). 

Definition  5.3:  A  mapping  is  a  sense-preserving  rigid  function  from  R-  onto  itself. 

Definition  5.4:  If  X  is  a  point  or  set  of  points  in  R^  and  M  is  a  mapping,  then 
image(Af,X)  =  Af(X). 

Definition  5.5:  A  general  velocity  is  the  derivative  of  a  function  from  the  reals  onto 
mappings. 

Definition  5.6:  An  object  is  a  primitive  entity. 

Definition  5.7:  For  any  object  0,  shape(O)  is  a  subset  of  R^. 

Definition  5.8:  A  scene  is  a  function  on  a  set  of  objects,  which  takes  each  object  into  an 
ordered  pair  of  a  mapping  and  a  general  velocity.  For  any  object  0  and  scene  S, 
mapping(5,0)  is  the  mapping  specified  by  S  for  object  O;  pnor-velocity(S,0)  is  the  general 
velocity  specified  by  S  for  0. 

Definition  5.9:  For  any  scene  S  and  object  O,  place(0,5)  =  mappingCS',  shape(O)),  the 
point  set  occupied  by  O  in  S. 

Definition  5.10:  An  object-state  is  an  ordered  triple  of  an  object,  a  mapping,  and  a 
movement.  For  any  scene  S  and  object  O, 
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state(0,5)  =   <  0,  mapping(5,0).  moveraent(5,0)  > 

Definition  5.11:  A  pseudo-object  is  an  ordered  pair  of  an  object  and  a  geometric  entity. 
We  write  Q  =  pseudo-obiect(0,A:X  ) 

Definition  5.12:  Let  Q   =   pseudo-object(O..YA')  be  a  pseudo-object.   Then  source(2)   = 
0,  and  shape(C)  =   PP     For  any  scene  5,  place(2,5)  =   image(mapping(5,0)./'P). 
Definition  5.13: 

A  chroriicle  is  a  function  from  an  interval  of  time  into  scenes.  For  any  chronicle  C. 
duration(C)  is  the  domain  of  C,  a  time  interval.  For  any  C  and  any  T  €  duration(C), 
scene(C,7')  =  C{T),  the  scene  of  chronicle  C  at  time  T. 

All  chronicles  are  subject  to  the  following  constraints: 

i.        The    duration   of   the    chronicle    is    closed    on    the    left.    Each    chronicle    has   a   starting 

moment, 
ii.       Any  two  scenes  in  the  range  of  the  chronicle  have  the  same  objects  in  their  domain. 

iii.  For  any  object  O  in  the  chronicle,  the  function  \(7")  mappingCscenefCD.O)  is 
continuous.  That  is,  objects  move  continuously  in  space. 

iv.  For  any  object  O  in  the  chronicle,  the  function  K{T)  prior-velocity(scene(C,r),0)  is 
continuous  from  the  left. 

V.  For  any  object  O  in  the  chronicle,  the  function  K(T)  prior-velocity{scene(C.7"),0)  is 
equal  to  the  derivative  from  the  left  of  the  function  k{T)  mapping(scene(C,r),0),  at  all 
T  after  the  start  of  the  chronicle.  That  is,  the  velocity  is  the  derivative  of  the  position. 

Note  that  chronicles  do  not  have  to  be  physically  possible.  We  will  use  the  predicate 
"phys-poss(C)"  to  distinguish  chronicles  that  obey  the  laws  of  physics. 

The  next  four  definitions,  relating  to  force,  are  not  needed  in  anything  that  follows. 
They  are  included  only  to  show  how  forces  can  be  handled  in  this  system. 

Definition  5.14:  A  force  is  a  vector  with  the  dimensions  of  force  (mass  times  distance 
over  time  squared). 

Definition  5.15:  A  force  distribution  is  a  function  f{XX)  from  a  measureable  point  set 
XX  to  a  force.  f(XX)  measures  the  total  force  exerted  across  the  volume  XX.  A  force 
distribution  must  satisfy  the  following  axiom: 

If  XXI  and.YXZ  are  disjoint,  then   f{XXl\jXX2)  =  f(XX])' f(XX2) . 

Definition  5.16:  A  force  account  is  a  function /(01,A', 02)  which  maps  two  objects  01. 
01  and  a  flag  A'  onto  a  force  distribution.  01  is  the  object  acted  on  by  the  force.  N  is  one  of 
the  four  constants,  "gravity,"  "normal,"  "friction."  or  "drag."  02  is  the  source  of  the  force, 
for  normal  and  friction  forces. 

Definition  5.17:  A  force  account  chronicle  is  a  function  from  instants  of  time  to  force 
accounts. 

Using  these  concepts,  we  can  state  the  axioms  necessary  to  do  force  analysis.  The  key 
predicate  here  is  '■compatible(FW,C) "  which  states  that  a  force  account  chronicle  FH  is 
compatible  with  a  chronicle  of  motion  C.  A  chronicle  is  physically  possible  only  if  there 
exists  some  compatible  force  account  chronicle.  The  basic  axioms  governing  compatability  are 

i.  Newton's  second  law  F  =  ma.  The  acceleration  of  an  object  in  the  chronicle  must  be 
proportional  to  the  total  force  exerted  on  it  in  the  force  account  chronicle. 

ii.  Newton's  third  law.  If  A'  is  a  friction  or  normal  force,  and  f"  is  a  force  account,  then 
F(0\,N  ,02)  =  F{02J^  ,0\). 

iii.  Particular  laws  governing  each  different  kind  of  force.  For  example,  normal  and  friction 
forces  act  only  at  the  surface  of  an  object. 
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6.    Axioms  for  Physical  Reasoning 

In  this  section  we  develop  a  first  order  language  and  axioms  adequate  to  solve  the  first 
"die  in  the  funnel"  example.  All  rules  will  be  stated,  first  in  English,  then  in  first-order 
logic.  We  use  a  ryped  logic,  indicating  the  type  of  a  variable  by  its  first  letter.  The  following 
types  are  used,  with  the  indicated  keys. 


Type 

First  letter 

Real-valued  quantity 

R 

Instant  of  time 

T 

Point  in  R^ 

X 

Vector 

V 

Vector  field 

p 

Mapping 

M 

Object 

O 

Pseudo-object 

Q 

Scene 

S 

Chronicle 

c 

Sets  of  entities  will  be  designated  by  a  double  initial  letters.  For  example,  variables 
denoting  intervals  of  time  will  begin  TT;  variables  denoting  sets  of  objects  will  begin  OO; 
etc.  On  occasion  it  will  be  useful  to  state  axioms  which  can  apply  either  to  objects  or  to 
pseudo-objects;  variables  of  this  kind  will  begin  with  Q. 

Capital  letters  are  used  for  variables;  lower-case  letters  for  predicates,  functions,  and 
constants.  Variables  not  explicitly  quantified  are  assumed  to  be  universally  quantified,  with 
the  entire  formula  as  scope. 

The  non-logical  terms  and  axioms  of  set  theory  and  real  analysis  are  assumed.  In  our 
English  definitions  of  functions  and  predicates,  we  will  assume  common  geometrical 
concepts.  Our  list  of  axioms  will  be  divided  into  geometrical  axioms,  axioms  of  motion,  and 
physical  axioms.  Definitions  of  predicates  ,ind  functions  will  be  interspersed  with  the 
statement  of  the  axioms  and  theorems  which  use  them.  The  predicate  calculus  formula 
following  the  English  statement  of  the  definition  may  be  equivalent  to  the  definition;  or  it 
may  be  some  weaker  statement  closely  associated  with  the  definition.  Subsequent  proofs 
using  the  definition  are  allowed  only  to  use  the  predicate  calculus  form;  though  subsequent 
axioms  may  be  based  on  the  full  definition.  Note  that  the  only  axioms  and  theorems  listed 
are  those  we  need  for  our  particular  example.  .Many  other  important  and  true  axioms  are 
omitted.  Theorems  in  this  section  are  statements  which  could  be  proven  from  the  axioms,  but 
which  are  envisioned  as  part  of  the  reasoner's  permanent  stock  of  rules.  No  proofs  of 
theorems  are  given  in  this  section. 

6.1.    Geometry 

6.1.1.    Vector  algebra. 

We  assume  the  use  of  addition  and  subtraction  of  vectors,  the  addition  of  vectors  to 
points  to  give  another  point,  the  subtraction  of  points  to  give  a  vector,  and  the  multiplication 
of  vectors  by  scalars.  We  also  assume  the  use  of  the  dot  product  function  (which  we  will 
write,  "V1V2")  mapping  two  vectors  into  a  real.    We  assume  all  the  usual  axioms  on  these. 

Definition  6.1.1.1:  The  function  "length(U) '  maps  a  vector  V  into  its  length. 
length(V)  •  length(V)  =   V  V 

length(V)  >  0 
.\xiom  6.1.1.2:  The  dot  product  of  rwo  vectors  is  at  most  the  product  of  their  lengths. 
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abs(VlV2)  <  length(Vl)     lengih(V2) 

Definition  6.1.1.3:  The  function  "distance(Xl,X2)'  maps  two  points  .VI  and  X2  onto  the 
distance  bet^^■een  them.  In  the  usual  way,  we  extend  the  distance  function  to  taice  closed  point 
sets  as  arguments,  meaning  the  minimum  distance  between  any  points  in  them. 

distance(Xl,X2)  =   length(.Yl- .Y2) 

distance(XA:i,XX2)  =   glb(distance(A:i,A'2)  |  A:1€A:X1,.Y2€XA:2) 

Axiom  6.1.1.4:  The  distance  between  two  points  is  0  iff  they  are  equal. 
distance(A:i,A:2)  =  0  »  X1  =  .V2 

Definition  6.1.1.5:  The  predicate  "paralIel(Vl,V2)"  holds  on  two  vectors  VI,  V2  if  they 
are  parallel. 

parallel(V'l,V2)  c=>  [V\  =  0  v  3^V2  =  /?  VI] 

Definition  6.1.1.6:  The  predicate  "beyond(A'l,.V2,V)"  (read  "XI  is  beyond  X2  in 
direction  V")  holds  for  points  XI,  X2  and  vector  V  if  (X1-X2)V  is  positive.  Applied  to 
point  sets,  in  the  form  "beyond(XXl,X.Y2.V)"  it  means  that  every  point  of  the  first  is  beyond 
every  point  of  the  second. 

beyond(Xl,X2,V)  o  (X1-X2)'V>0 

beyond(XXl,XX2,V)  o  [  "i xi-xx]  .v:..v,y:  beyond(Xl,X2,V)  ] 

Definition  6.1.1.7:  The  predicate  "weak-beyond(Xl,.X2,V)"  fread  "XI  is  weakly  beyond 
X2  in  the  direction  of  V")  holds  for  points  XI,  X2  and  vector  V  if  (X1-X2)V  is  non- 
negative.  Applied  to  point  sets,  in  the  form  "weak-beyond(XXl,XX2,V)"  it  means  that  every 
point  of  the  first  is  weakly  beyond  every  point  of  the  second. 

weak-beyond(Xl,X2,V)  o  not(  beyond(X2,Xl,V)) 
weak-beyond(XXl,.VX2,V)  <=>  [  '^ x\-.x.\\.X2iXX2  weak-beyond(Xl.X2,V)  ] 

6.1.2.    Basic  Concepts 

Definition  6.1.2.1:  The  function  "make-line(X,V)"  is  the  line  through  point  X  parallel  to 
vector  V. 

XHmake-line{X,V)  o  parallel((Xl-X),V)  ] 

Definition  6.1.2.2:  The  predicate  "line(XX)"  holds  if  the  point  set  XX  is  a  line. 
line(XX)o3vi   V^O  '  XX  =  make-Hne(X,V) 

Definition  6.1.2.3:  The  predicate  "parallel-lines(XXLl,XXL2)"  holds  if  XXLl  and  XXZ.2 
are  parallel  lines. 

parallel-lines(XX/:i,XXL2)  o 

3v,.v:.i  XXZ.l  =  makeTine(Xl,V)  a  XXL2  =  make-line(X2,V) 

Theorem  6.1.2.4:  (Parallel  Postulate)  Given  a  point  X  and  a  line  XXL,  there  is  a  line 
parallel  to  XXL  through  X. 

3v.Yi:  paralIel-lines(XXZ.2,XXZ.)  a  X^XXLl 
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Definition  6.1.2.5:  The  function  'line-seglA' 1  ,X2)"  maps  two  points  XI.  X2  onto  the  line 
segment  connecting  them. 

X€line-seg(Xl,A'2)  ^  3,  0</?<l   •.  X  =  Xl~  R(X2-Xl) 

line-seg(A'l.A'2)  C  make-line(Al.A2-Xl) 

Definition  6.1.2.6:  The  predicate  "single-point(XA)"  holds  if  point  set  XX  consists  of  a 
single  point. 

single-point(AA)  <=>  3s^XX  =  {X\ 

Definition  6.1.2.7:  The  function  "interior(AA)"  maps  a  point  set  XX  onto  its  interior. 
"Interior"  is  interpreted  in  the  sense  of  a  manifold;  a  point  XdXX  is  in  the  interior  of  XX  if 
the  neighborhood  of  X  in  XX  is  homeomorphic  to  ^"  for  some  k.  Thus,  lines,  surfaces,  and 
solids  all  have  their  own  natural  sense  of  interior. 

interior(AX)  C  XX 

Definition  6.1.2.8:  The  function  "boundary(AX)"  maps  a  point  set  XX  onto  its  boundary. 

The  boundary  of  XX  is  equal  to  its  closure  minus  its  interior  (as  defined  above). 

boundary(A'A)  H   mterior(AX)  =  0 

Definition  6.1.2.9:  The  predicate  "body(XA)'  means  that  XX  is  the  closure  of  a 
connected,  non-null,  three  dimensional  manifold. 

Definition  6.1.2.10:  The  predicate  "surface(A-X)"  means  that  AX  is  the  closure  of  a 
connected,  non-null,  two  dimensional  manifold. 

Axiom  6.1.2.11:  A  body  or  a  surface  is  the  union  of  its  boundary  and  its  interior. 
[  body(XX)  \.  surface(XX)  ]  =>  XX  =  boundary(XX)  |J  interior(XX) 

Definition  6.1.2.12:  The  predicate  "abut(XXl,XX2,XX/i )"  is  true  if  XXI  and  XX2 
intersect  only  in  their  boundaries,  and  XX/\  is  their  intersection. 

abut(XXl,XX2,XX/^)  o  XXA   =  XX\P\XX2  =   boundary(XXl)  H   boundary(XX2)  *  <Z) 

Definition  6.1.2.13:  The  predicate    ■overlap{XXl.XX2)"  holds  of  point  sets  XXI,  XX2  if 

their  interiors  intersect. 

overlap(XXl,XX2)  <=>  interior(XXl)  H   interior(XX2)  i^  0. 

Theorem  6.1.2.14:  If  two  point  sets  abut,  then  they  do  not  overlap, 
not  [  abut(XXl,XX2,XX/\)  a  overlap(XX  1,XX2)  ] 

Axiom  6.1.2.15:  Any  two  bodies  either  overlap,  abut,  or  are  disjoint. 

[  body(XXl)  A  body(XX2)  ]  => 

[  overlap(XXl,XX2)  v  3^^^  abut(XXl.XX2,XX/l )  v  disjoint(XXl,XX2)  ] 

Definition  6.1.2.16:  The  predicate  "convex(XX)"  means  that  the  point  set  XX  is  convex 
in  the  usual  sense:  any  two  points  in  XX  can  be  connected  by  a  line  segment  lying  in  XX 

convex(XX)  ^  '^ xxxz.xx  line-seg(Xl,X2)CXX 

Definition  6.1.2.17:  The  function  "convex-hull(XX)"  maps  a  point  set  XX  to  its  convex 
hull. 

convex(convex-hull(XX))  a  XX  C  convex-hull(XX) 
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Axiom  6.1.2.18:  A  convex  point  set  is  its  own  convex  hull. 
convex(A'A')  =>  XX  =  convex-hull(A'.V ) 

Definition  6.1.2.19:  The  predicate  ■strongly-convex(XX ) '  means  that  any  two  points  m 
XX  can  be  connected  by  a  line  segment  lying  in  the  interior  of  XX. 

strongly-convex{XX)  o 

'''yia-ax  (line-seg(.X:i,X2)-.Yl-.V2)  C  interior(A'A') 

Theorem  6.1.2.20:  A  strongly  convex  shape  is  convex. 
strongly-convex(AA)  =>  convex(AA) 

Theorem  6.1.2.21:  If  a  strongly-convex  shape  abuts  a  convex  shape,  they  abut  in  a 
single  point. 

[  strongly-convex(AA'l)  a  convex(A:X2)   ■-  abut(A'A  l,AA2,AAf )  ]  =. 
single-point(AXF) 

Definition  6.1.2.22:  The  function  "diameter(AA )"  maps  a  point  set  XX  onto  the 
maximum  distance  between  any  two  points. 

diameter(XX)  =  lub(distance(Al,X2)  \Xl,X2£XX) 

Theorem  6.1.2.23:  If  point  set  XXI  intersects  both  point  set  XX2  and  point  set  XX3, 
then  the  distance  from  XX2  to  XX3  is  less  than  the  diameter  of  XXI. 

[  XXin'^^2^0  A  XXinA:X3^0  ]  i>  distance(.VX2.XX3)  <  diameter(XXl) 

Definition  6.1.2.24:  The  predicate  "congruent(XXl,XX2)"  holds  for  two  point  sets  XXI. 
XX2  which  are  congruent  without  reflection;  that  is,  there  is  a  rigid  mapping  from  XXI  to 
XX2. 

congruent(XXl,XX2)  «  3,,  XX2=  image(M,XXl) 

Definition  6.1.2.25:  The  predicate  "bounded(XX,V')"  means  that  point  set  XX  is  bounded 
in  direction  V. 

bounded(XX,U)  ^  3^|V^^.y  beyond(Xl,X,V) 

Definition  6.1.2.26:  The  predicate  "extreme(X.XX,V')"  means  that  X  is  an  extreme  point 
of  XX  in  direction  V. 

extreme(X,XX,V)  <=>  [X€XX  a  V  \.-,  .^^  weak-beyond(X,X2,V)  ] 

Axiom  6.1.2.27:  If  a  body  is  bounded  in  direction  V  then  it  has  an  extreme  point  in 
direction  V. 

[  body(XX)  A  bounded(XX,V)  ]  ==>  3  ^  extreme(X,XX,V) 

Axiom  6.1.2.28:  An  extreme  point  of  a  body  is  on  its  boundary. 
[  extreme (X, .XX, V)  a  body (XX)  ]  =>  Xe boundary (XX) 

6.1.3.    Planar  Geometry 

Definition  6.1.3.1:  The  function  "make-plane(X,V)"  gives  the  plane  through  X  normal  to 
V. 

XHmaiie-plane(X,V)  o  (Xl-X)  V  =  0 
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Definition  6.1.3.2:  Ihe  predicate  "plane(.VA')"  holds  if  XX  is  a  plane. 
plane(A'X)  ■==>  3y.    V^O  ■  XA' =  make-plane{X.V) 

Definition  6.1.3.3:  The  predicate  "planar(A'A')"  holds  if  point  set  A' A  lies  in  a  plane. 
planarC.YA)  «  3  ^^.^.^  plane(AAT)   ■  XX  C  XXP 

.4xiom  6.1.3.4:  A  point  and  a  line  are  coplanar. 
Iine{AA-L)  =>  planarl.VAZ.  (J  A) 

Definition  6.1.3.5:  The  predicate  cross-sectionCAACA.YB  ,.VA/')"  means  that  AAC  is  the 
closure  of  a  connected  component  of  the  intersection  of  the  interior  of  body  AA'O  with  plane 
XXP.  (This  convoluted  definition  is  necessary  to  handle  cases  where  plane  cuts  a  body  in 
several  surfaces,  where  a  plane  is  tangent  to  a  body,  or,  worse,  where  it  is  tangent  for  a 
while  and  then  enters  the  body.) 

cross-section(AAC.AAfi,AAP)  => 

[bodyCVAfi)   ■    plane(AAP)   '   surface(AAC)   -  XXC CXXP  {~]XXB  ] 

.\.\iom  6.1.3.6:  .■\ny  interior  point  of  body  AAfi  in  plane  XXP  is  in  some  cross  section  of 
XXB  by  XXP . 

[  body(AAfi)   ■   planeCAAP)  /  A€  interior(AAfl)  H^^^  1  => 
3  xvc  ^  i  'VAC  •  cross-section(AAC  .XXB  .XXP ) 

Axiom  6.1.3.7:  The  boundary  of  a  cross  section  lies  in  the  boundary  of  the  solid. 
cross-sectionfXACAAfi.AXP)  =>  boundary (.V.YC)  C  boundary(AAB) 

Theorem  6.1.3.8:  If  two  bodies  do  not  overlap,  then  no  cross  sections  by  any  plane 
overlap. 

[  not(overlap(AA01,AA02))   ■  cross-section(A:ACl,.'!:A01,A:X/') 

/  cross-section(AXC2,AX02.AAP)  ]  =. 
not(overlap(AAC  1  .AAC2)) 

Axiom  6.1.3.9:  .Any  two  co-planar  surfaces  either  overlap,  abut,  or  have  null 
intersection. 

[  planar(AAlUAA2)  /   surface(AAl)  '   surface(AA2)  ]  ^ 

[  overlap(AAl,AA2)  v  3  ^.^.^  abut(AAl,AA2,.\:AA  )  v  disjoint(AAl,X.Y2)  ] 

Definition  6.1.3.10:  The  function  "top(AA,K)"  maps  a  closed,  bounded,  point  set  AA 
and  a  vector  V  into  the  plane  perpendicular  to  V^  passing  through  the  extreme  point(s)  of  AA 
in  the  V  direction.  The  function  "bottom(AA,V)"  is  the  same  in  the  -V  direction. 

extreme(A,A:A',V)  =>  topiXX ,V)  =  make-plane(A,V) 
bottom(AA,V)  =  top(AA,-V) 

Definition  6.1.3.11:  The  predicate  "perpendicular(V',AX/')"  means  that  vector  V  is 
perpendicular  to  plane  AAP. 

perpendicular(V,AA/')  <s>  [  plane{XXP)  a  V  v,  .v:^,v.vp  V'(A1-A2)  =  0  ] 

Definition  6.1.3.12:  The  predicate  "pparallel(V,.YX/')"  means  that  vector  V  is  parallel  to 
some  vector  in  plane  XXP. 
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pparalleUVXYP)  <=>  3,^.  x:,xxp  ^^*^-  '  parallel(V..Yl-X2) 

6.1.4.  Curve  characteristics 

Definition  6.1.4.1:  The  predicate  "sidc(XXS ,XXP)"  holds  for  point  sets  XXS  and  XXP  if 
XXP  is  a  planar  surface  and  XXS  is  a  connected  subset  of  the  boundary  of  XXP . 

sidc{XXS,XXP)  =>  planarlXXP)  a  surface(XXP)  a  A:a:5C boundary (XX/') 

Definition  6.1.4.2:  The  predicate  "p-smooth(XA',X)"  is  true  if  XX  is  a  planar  surface,  X 
is  in  the  boundary  of  XX,  and  the  boundary  of  XX  is  smooth  at  X. 

p-smooth(XX,X)  =>  planar(XX)  a  surface(XX)  a  X  €  boundary(XX) 

Definition  6.1.4.3:  The  function  "p-norm(XX,X)"  gives  the  unit  normal  to  the  boundary 
of  the  planar  surface  XX  at  the  boundary  point  X  directed  out  of  XX.  p-norm(XX,X)  is 
defined  only  if  the  boundary  of  XX  is  smooth  at  X. 

length(p-norm(XX,X))  =   1 

Definition  6.1.4.4:  The  predicate  "d-sidciPC ,XX)"  (read  "PC  is  a  directed  side  of  XP") 
is  true  if  XX  is  a  planar  surface,  and  PC  is  the  vector  field  of  the  normal  to  XX  along  some 
side  of  XX. 

d-s\dc(PC,XX)  ■» 

[  side(base(PC),XX)  ■  V  ^.  Xihase(PC)  =>  value(PC.X)  =  p-norm(XX,X)  ] 

Definition  6.1.4.5:  The  predicate  "convex-side(P  ,XX)"  holds  for  a  vector  field  P  if  P  is 
the  directed  side  of  XX  and  also  of  some  convex  shape. 

convex-sidelP.XX)  <=> 

[  d-side(P,XX)  a  3^.^.,  convex(XXl)  a  d-sidefP  ,XX1)  ] 

Axiom  6.1.4.6:  If  a  strongly  convex  surface  abuts  another  co-planar  surface  along  a 
convex  side  of  the  latter,  then  they  abut  in  a  single  point. 

[  strongly-convex(XXl)  a  abut(XXl,XX2,XX-4)  a  planar(XXHJXX2)  a 

surface(XXl)  a  surface(XX2)  a  XXACbase(P)  a  convex-sidetP  ,XX2)  ]  => 
single-point(XX4 ) 

Definition  6.1.4.7:  The  function  "inner-dsidefXXf  ,XXZ.)"  takes  as  arguments  a  planar 
object  XXP  and  a  disjoint  line  XXL  in  the  same  plane.  Its  value  is  defined  as  follows:  Drop 
the  perpendiculars  from  the  top  and  bottom  of  XXP  to  XXL.  These  divide  the  boundary  of 
XXP  into  at  least  rwo  parts.    The  part  which  lies  closest  to  XXL  is  the  inner  side.  (Figure  24) 

[  planarCXXP^jXXZ.)  a  XXLf]XXP  =  0  a  line(XXZ.)  ]  => 
d-side(inner-dside(XXP,XXL),XX/') 

Axiom  6.1.4.8:  If  an  inner  side  of  a  shape  with  respect  to  a  line  is  convex  and  smooth, 
then  the  normal  to  the  side  is  never  parallel  to  the  line. 

[  /'5  =  inner-dside(XX/',make-line(X0,V))  a  convex-side(P5,XXP)  a  Xibase{PS)  ]  => 
[  not(parallel(value(P5,X),V))  a  not(parallel(value(/'5,X),- V))  ] 

6.1.5.  Surface  Characteristics 

Definition  6.1.5.1:  The  predicate  "face(XXf  ,XXO) "  holds  for  point-sets  XXF  and  XXO  if 
XXO  is  a  body,  and  XXF  is  a  connected  subset  of  boundary(XXO). 
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iaceiXXF .XXO)  ^  [  XXf  CboundaryfX.VO)   -  XXF^CS  ] 

Definition  6.1.5.2:  The  predicate  "smooth(,YX,A' )"  is  true  if  XX  is  a  body,  A'  is  in  the 
boundary  of  XX ,  and  the  boundary  of  XX  is  smooth  at  X . 

smooth(A:A:,A:)  =>  body(XX)  '  X  €  boundary (XA) 

Definition  6.1.5.3:  The  predicate  "ordinary(XX)"  is  true  if  XX  is  a  body  which  is 
everywhere  smooth. 

ordinary(XX)  o 

[  body(XX)  A  V\.  [  X«boundary(XX)  =>  smooth(XX,X)  ]  ] 

Definition  6.1.5.4:  The  function  "surf-norm(XA',X)"  gives  the  unit  surface  normal  to  a 
body  XX  at  a  point  X  on  the  boundary  of  XX.  surf-norm(XX,X)  is  defined  only  if  XX  is 
smooth  at  X .    surf-norm(XX  ,X )  is  directed  out  of  XX . 

Axiom  6.1.5.5:  Smoothness  and  the  value  of  the  surface  normal  are  local  properties  of 
the  boundary.  Specifically,  if  two  bodies  share  part  of  their  boundary,  then,  at  any  interior 
point  of  the  overlap,  one  is  smooth  iff  the  other  is  smooth,  and  their  surface  normals  are 
either  parallel  or  anti-parallel. 

[body(XXl)  ^  body(XX2)  a  XXA  C  boundary(XXl)  f]  boundary(XX2)  a 

and  X€interior(XX/\)  a  smooth(XXl,X)  ]  => 
[  smooth(XX2,X)  a 

[  surf-norm(XXl,X)  =  surf-norm(XX2,X)   y  surf-norm(XXl,X)  =   -surf-norm(XX2.X) 

Axiom  6.1.5.6:  If  two  bodies  have  a  common  smooth  boundary  point,  and  the  surface 
normals  at  the  point  are  not  anti-parallel,  then  the  bodies  must  overlap. 

[  body(XXl)  A  body(XX2)  .    X€  boundary(XXl)  f]  boundary(XX2)  a 

smooth(XXl,X)  ^  smooth(XX2.X)  ]  => 
[  surf-norm(XXl,X)  =   -  surf-norm(XX2,X)  .  overlap(XXl,XX2)  ] 

Theorem  6.1.5.7:  If  bodies  XXI  and  XX2  both  abut  XX3  at  the  same  smooth  boundary 
point  of  all  three  surfaces,  then  XXI  overlaps  with  A'X2. 

[  abut(XXl,XX3,XXA)  a  abut(XA2,XX3,XXB)  '  X^XXAf)XXB  a 

smoothfXXl.X)  '  sraooth(XX2,X)  /  smooth(XX3,X)  ]  =. 
overlap(XXl.XX2) 

Axiom  6.1.5.8:  If  the  extreme  point  of  a  body  m  direction  V  is  a  smooth  point  of  its 
boundary,  then  the  surface  normal  there  is  parallel  to  V . 

[  smooth(XX,X)  A  extreme(X,XX,V))  ]  =>  parallel(surf-norm(XX,X),V) 

Axiom  6.1.5.9:  If  a  convex  shape  is  smooth  at  a  boundary  point  X,  then  X  is  an  extreme 
point  of  the  shape  in  the  direction  of  the  surface  normal  at  X. 

[convex(XX)  a  smooth(XX,X)]  =>  extreme(X,XX,surf-norm(XX,X)) 

Axiom  6.1.5.10:  Given  a  smooth  boundary  point  on  a  body,  construct  a  plane  through 
the  point  containing  the  surface  normal  at  the  point.  Then  (i)  there  is  a  cross  section  of  the 
body  by  the  plane  containing  the  point;  (ii)  the  cross  section  is  smooth  at  the  point;  and  (hi) 
the  normal  to  the  cross  section  is  equal  to  the  surface  normal  at  the  point. 

[  body(XX)  A  smooth(XX,X)  a  pparallel(surf-norm(XX,X),XX/')  ]  => 
[  3yvc  cross-section(XXC,XX,XXP)  a  X€XXC  a 
p-smooth(XXC,X)  a  p-nonn(XXC,X)  =  surf-norm(XX,X")  ] 


Definition  6.1.5.11:  The  function  "const-face (.VA'.V)'  maps  a  planar  shape  XX.  and  a 
vector  V  into  the  vector  field  equal  to  V  at  every  point  of  XX . 

base(const-face(A:A:,V))  =  XX 

X£XX  =>  vaiuelconst-facelXX.V),^)   =   V 

Definition  6.1.5.12:  The  function  "dboundary(.VX)"  (read,  "the  directed  boundary  of 
XX")  maps  an  ordinary  point  set  XX  onto  the  vector  field  which  points  outwards  from  XX  at 
each  boundary  point. 

ordinary(XX)  => 

[  base(dboundary(XX))  =  boundaryCXX)  ^ 
Vv  X€boundary(XX)  =>  value(dboundary(XX),X)  =   surf-norm(XX,X)  ] 

Definition  6.1.5.13:  The  predicate  "diaceiP ,XX)"  holds  for  a  vector  field  P  and  a  point 
set  XX  if  hasc{P)  is  a  face  of  XX,  XX  is  smooth  throughout  P,  and,  for  all  X^P ,  value(/'.X) 
=  surf-norm{XX,X). 

dface(F,XX)  « 
[  face(base(/'P),X)  a 
Vv[Xebase(P)  =>  [  smooth(XX,X)  a  value(P,X)  =  surf-norm(XX.X)]]] 

Definition  6.1.5.14:  The  predicate  "in-front(X,P )"  holds  for  a  point  X  and  a  directed 
face  P  if,  for  some  point  in  X 1  in  P ,  X  is  beyond  X 1  in  the  direction  of  P  at  X 1 . 

in-front(X,/')  <=>  3^,  XKbase(P)  a  beyond(X,Xl,value(/',Xl)) 

Definition  6.1.5.15:  The  predicate  "weaji-in-frontfX.P)"  holds  for  a  point  X  and  a 
directed  face  P  if,  for  some  point  in  XI  in  P,  A  is  weakly  beyond  XI  in  the  direction  of  P  at 
XI. 

weak-in-front(X,P)  <=■  3^.,  XKbase(P)  a  weak-beyond(X,X],value(P,Xl)) 

Definition  6.1.5.16:  The  predicate  "f-abut(P  1,P2)"  holds  of  two  vector  fields  PI,  P2  if 
PI  and  P2  are  directed  faces  of  two  shapes  .XXI  and  XX2,  and  XXI  and  XX2  abut  in  the 
intersection  of  PI,  P2. 

f-abut(Pl,P2)  o 

^xwxxz.xx.-^  '^^^^^(P'^-^^'^)   '  dface(P2,.YX2)   ■    abut(.XXl,XX2,XXA )  a 
XXA  C  base(Pl)  p)  base(P2) 

Axiom  6.1.5.17:  If  the  boundary  of  a  convex  shape  f-abuts  a  vector  field,  then  the  shape 
lies  weakly  in  front  of  the  field. 

[  f-abut(dboundary(XX),P)  a  convex(XX)  ]  ^  weak-in-front(XX,P) 

6.1.6.    Topology 

Definition  6.1.6.1:  The  predicate  "box(XX)"  holds  if  XX  is  an  ordinary  point  set  and  the 
complement  of  XX  has  exactly  two  connected  components. 

box(XX)  =>  ordinary (XX) 

Definition  6.1.6.2:  The  functions  "insideC.VX)"  and  "outside(XX)"  are  defined  if  X.X  is  a 
box,  and  denote  the  inside  and  outside  of  the  box.  Technically,  inside(XX)  is  the  closure  of 
the  finite  component  of  the  complement  of  XX;  outside(XX)  is  the  closure  of  the  infinite 
component  of  the  complement  of  XX. 
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Axiom  6.1.6.3:  The  inside  and  outside  of  a  box  abut  the  box. 

boxCXX)  => 

[  3vx4XVB  abut(inside(A'X).A:X.A:X>l)    -  abut(outside(.V.Y  ).A'A',.VA'fl )  ) 

Definition  6.1.6.4:  The  functions  "inner-boundary(.VX)"  and  "outer-boundarylA'A)"  are 
defined  if  XX  is  a  box.  They  denote  the  inner  and  outer  boundaries  of  XX . 

box(XX)  => 

[  inner-boundaryCA-A)  =  inside(.VA)   H  ^'^  ■'  outer-boundary(AA)  -  outside(A'A)  H  ^^ 

Definition  6.1.6.5:  The  predicate  'onficef.X'AO.AAS )"  holds  if  XXO  is  a  body,  XXB  is  a 
box,  and  XXO  is  a  subset  of  XXB  which  goes  from  the  inside  to  the  outside. 

orifice(AAO,AXS)  o 

[  box{XAfi)   'X  body(AAO)  a  XXO  C  XXB  ^ 
face(A:A:oninside(AAfi),  inside(AA:5))   '   face(AAO  H  outside(AAfl  ),  outside(AAB))  ] 

Definition  6.1.6.6:  The  predicate  "lubeiXXT ,Pl.P2y  (read  "AAr  is  a  tube  with  caps  PI 
and  P2")  if  PI  and  P2  are  directed  faces  which  "close  off"  AAT  into  a  box.  Formally,  PI  and 
P2  are  the  directed  faces  of  two  disjoint  solids  XXCl  and  XXC2  which,  together  with  .V.^T", 
form  a  box  XXB  in  such  a  way  that  each  of  AA'Cl,  XXC2,  and  AAT"  is  an  orifice  of  the  box. 
(Figure  25) 

tube{XXT,P\,P2)  o 

3 y^c ,  VVC2 .XXB  [  ^^B  =  -^-^^  U  ^^'^  1  U  •^'•^'<^ -   '  box(AAfi )  .■ 

orifice(A'ACl.AXS)  a  orifice(AAC2,A'A:fl )  a  orifice(AXr,AXfl)  a 

base(Pl)  =  boundary(AACl)  p)  outer-boundary(AAS)  a 

base(P2)  =  boundary(AAC2)  f]  outer-boundary(AXfl)  a 

df  ace  (PI, XXCl)  a  dface(P2,XXC2)   •  disioint(XXCl,XXC2)  a 

not(overlap(XXCl,XXr))  a  not(overlap(XXC2,XXr))  ] 

Definition  6.1.6.7:  The  function  'tube-inside(XXr,Pl,P2) '  applied  to  a  tube  XXJ  with 
caps  PI,  P2  gives  the  (closed)  volume  inside  the  tube  defined  by  XAT,  PI.  P2. 

PI\JP2  C  dboundary(tube-inside(XAr.Pl.P2)) 

boundary(tube-inside(XX7,Pl.P2))  -  base(Pl)  -  base(P2)  C  boundary(XX7") 

Definition  6.1.6.8:  The  function  "inner-surface(.X'.\T,Pl,P2) '  maps  a  tube  XAT  with  caps 
P  1 ,  P2  into  the  part  of  the  boundary  of  AAr  on  the  inside. 

tube(XX7",Pl,P2)  ^ 

inner-surface(XXr,Pl,P2)  =  boundary(XXr)  f]  boundary(tube-inside(XX7",Pl,P2)) 

Axiom  6.1.6.9:  If  XXT  is  a  tube  with  caps  PI,  P2,  then  the  inner  surface  of  XXr  abuts 
PI  and  P2  along  their  entire  boundaries.  These  rwo  boundaries  together  make  up  the  entire 
boundary  of  the  inner  surface. 

[  tube(XXr,Pl,P2)  A  XXI  =  inner-surface(XX7",Pl,P2)  ]  => 

[  abut(XX/,base(Pl),boundary{base(Pl))  a  abut(XX/,base(P2),boundary(base(P2))  a 
boundary(XX/)  =  boundary(base(Pl))  IJ  boundary(base(P2))  ] 

Definition  6.1.6.10:  The  predicate  •p-tube(XXr,Pl.P2,V) '  holds  if  XXT"  is  a  tube  with 
caps  PI,  P2,  where  PI  and  P2  lie  in  the  top/bottom  of  XXT"  in  direction  V,  and  have  vectors 
equal  to  V/-  V. 
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p-tube(XXr,/'l,/'2.V)  .-=• 

i  iube(XXT.Pl.P2)       base(/'l)  C  top(A:AT,V)   ■   base(P2)  C  bottom(.YA'r.V')  / 
P\  =  const-face(base(/'l).V  )      P2  =  const-face(base(/'2),  -  V)  | 

Axiom  6.1.6.11:  The  inside  of  a  p-tube  lies  entirely  to  one  side  of  either  cap. 

p-iuhe(XXT,Pl,P2,V)  => 

[  weak-beyond(tube-inside(.VX7',Pl,P2),base(/'l),- V)  a 
weak-beyond(tube-inside(.YX7  .Pl.P2),basc{P2).V)  ] 

6.1.7.    Special  Shapes 

Definition  6.1.7.1:  The  predicate  "circle(.Y.Y )"  means  that  XX  is  a  circle. 
Axiom  6.1.7.2:  A  circle  is  strongly  convex  and  planar. 

circle(XX)  =>  planar(.Y.Y)   \  strongly-convex(A'X) 

Definition  6.1.7.3:  The  predicate  "sphere(A'X) "  asserts  that  XX  is  a  sphere. 
Axiom  6.1.7.4:  A  sphere  is  strongly  convex. 
sphere(XX)  =>  strongly-convex(XX) 

Definition  6.1.7.5:  The  function  "radius(XX)"  is  defined  if  XX  is  a  sphere,  and  gives  its 
radius.  The  function  "center(XX)"  is  defined  if  XX  is  a  sphere,  and  gives  its  center. 

Axiom  6.1.7.6:  A  sphere  is  the  locus  of  all  points  within  the  radius  of  the  center. 
sphere(XX)  =>  [  Vv  X€XX  «  distance(X,center(XX))  <  radius(XX)  ] 

Theorem  6.1.7.7:  The  radius  of  a  sphere  is  half  its  diameter. 
sphere(XX)  =.  radius(XX)  =  1/2  •  diameter(XX) 

Axiom  6.1.7.8:  An  intersection  of  a  plane  and  a  circle  is  either  a  circular  cross  section 
or  a  single  boundary  point. 

[  sphere(XX5)  a  plane(XXP)  a  XXC  =  XXS  OXXP  i^0  ]  => 
[  [  single-point(,YXC)   '   XXC  C  boundary (XX5  )  ]   . 
[  circle(.VXC)      cross-section(XXC,XX5,XX/' )  |  ] 

Axiom  6.1.7.9:  The  center  of  a  sphere  is  an  interior  point. 
sphere(XX)  =>  center(XX)  €  interior(XX) 

Definition  6.1.7.10:  The  function  "solid-of-revolutionlX.Xf  .XXL)"  is  defined  if  XXL  is  a 

line;  it  is  the  solid  swept  out  rotating  XXF  around  XXL. 

Axiom  6.1.7.11:  A  sphere  is  a  solid  of  revolution  around  any  line  through  its  center. 
sphere(XX)  a  Vi=0  =>  3 v^c  XX  =  solid-of-revolution(XXC,make-line(center(XX),V)) 

Axiom  6.1.7.12:  If  two  solids  of  revolution  about  distinct  parallel  lines  abut,  then  the 
points  of  abutment  lie  in  the  plane  defined  by  the  center  lines. 

[  XXl  =  solid-af-revolution(XXf  l,XXZ.l)  a  XX2=  solid-of-revolution(XXF2.XXL2)  a 

XXLli^XXLl  A  parallel-lines(XXLl,XXZ.2)  a  abut(XXl,XX2,XX/4 )  ]  => 
planar(XXA  (JXXZ.1  (JXXi.2) 

Definition  6.1.7.13:  The  predicate   's-tube{XX)"  (read  "XX  is  an  s-tube")  holds  if  XX  is  a 
body  and  XX  is  a  solid  of  revolution  of  XXF  around  XXL  where  XXF  does  not  intersect  XXL. 
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s-tube(XX)  «. 

[  bodyCA-X)      BxxFX-c  ^irieiXXL)  ■  dis]oint(A:A:/",XXZ.)  a 
XX  =  solid-of-revolution(A:;f/",A:A:L)  ) 

Definition  6.1.7.14:  The  function  "center-line(XX)"  is  defined  if  XX  is  an  s-tube,  and 
returns  the  unique  line  around  which  XX  is  formed. 

XXZ,  =  center-line(XX)  <^  [  s-tube(XX)  and  3  ^^-^XX  =  solid-of-revoiution(XXF,XXL)  ] 

Definition  6.1.7.15:  The  function  "s-tube-top(XXF,V)"  maps  an  s-tube  XXF  and  a  vector 
V  along  the  center-line  of  V  into  a  cap  of  XAT  which  in  the  top  plane  of  XXF  in  direction  V. 
The  function  "s-tube-bot(XXf  ,V)"  is  analogous.  The  function  "s-tube-insidetXXF)"  gives  the 
tube-inside  of  XXf  with  respect  to  the  top  and  the  bottom. 

[  s-tube(XXf )  A  make-line(XO,V)  =  center-line(XXF)  a 

/'7  =  s-tube-top(X.YF,V)    a  PS  =  s-tube-bot(XXf  ,V)  ]  => 
[  base(Pr)  C  top{XXF,V)   \  Pr=  ccnst-face(base(/'),V)  a 

p-tube{XXF.PT.PB,V)  \  s-tube-inside(XXf )  =  tube-inside(XXf  .Pr.Pfl )  ] 

s-tube-bot(XXf  ,V)  =  s-tube-top(XXf ,- V) 

."Vxiom  6.1.7.16:  If  .VAT  is  an  s-tube  generated  from  a  planar  shape  XXF  around  a  line  in 
the  plane  of  XXF ,  then  the  intersection  of  XXT  by  any  plane  through  the  center  line  gives  two 
cross-sections.  Both  cross  sections  are  congruent  to  XXF  through  mappings  that  leave  XXL 
unmoved.  The  distance  ber\veen  them  is  equal  to  twice  the  radius  of  the  s-tube.* 

[  s-tube(XX7)   ■  XXT  =  soiid-of-revolutionfXXf  ,XXZ.)  a 

planar(XXFlJXXZ.)  a  planeiX.VP)       XXLCXXP  ]  => 
3.Y.YCI  v,YC2  wi  './:  [cross-section(XXCl..'CX7",XX/')  a  cross-section(XXC2,XXr,XX/')  a 

XXrnxX/'  =  XXCHjXXC2  /,  distance(XXCl,XXC2)  =  2 •  distance (XXf.XXZ.)  a 

XXCl  =  image(Ml,.YXf )   ■  XXL  =  image(Ml,XXL)  a 

XXC2  =  image(Af2,XXf )  ^  XXL  =  image (.V/ 2, XXL)  ] 

Axiom  6.1.7.17:  If  XXT  is  an  s-tube  generated  from  a  planar  shape  XXF,  then  any  radial 
cross  section  of  XX7  intersects  the  inner  boundary  of  XXr  in  a  curve  congruent  to  the  inner 
side  of  XXF. 

[  s-tube(XX7)  ,'  XX7  =  solid-of-revolution(XXF,XXL)  '  xXL  =  make-line(XO,V)  a 
planar(XXF(jXXL)  a  plane(XXF)      XXLCXXP   ^ 
XXIB  =  inner-boundary(XXr,s-tube-top(XXr,V),s-tube-bot(XX7',V))  a 
cross-section(XXC..VXr,XXF)  ]  => 

congruentiXXCpl-'i-'-'s^/fl.  base(inner-dside(XXF,XXL))) 

Axiom  6.1.7.18:  A  surface  normal  to  a  solid  of  revolution  lies  in  the  radial  plane. 

[  XX5'  =  solid-of-revolution(XXF,XXL)  \  smooth(XXS,X)  a 

plane(XXF)  a  X€XXF  a  XXLCXXP  ]  => 
pparallel(surf-norm(XXS,X),XXF) 


"This.  IS  an  axioml"  you  ask?  It  all  depends  where  you  want  to  invest  work.  If  we  had  bothered  to 
assert  something  about  the  solid-of-revolution"  function,  this  statement  could  now  be  proven.  But  it's 
easier  this  way.  As  Bertrand  Russell  said,  this  has  Uie  advantages  of  theft  over  honest  toil. 
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6.1.8.  Invariance  under  mapping 

All  the  functions  and  predicates  so  far  have  been  defined  without  reference  to  any  fixed 
coordinate  system.  Therefore,  they  are  all  invariant  under  rigid  mappings.  I  hat  is,  for  any 
function /(A  1,A2  ■  •  ■  Ak),  we  have  the  axiom 

image{MJ{Al,A2  ■  ■  ■  Ak))  =  /(image(A/ .-4  1).  image(Af  .,42)    ■         imagc{\f  .Ak]) 

For  each  predicate  p(Al,A2,  ■  ■  ■  Ak),  we  have  the  axiom 

p  (A  1,A2,  ■■■  Ak)  <:>  piimagsiM , A  I),  image(M ,A2)         •    image(A/ ,AA)) 

We  will  not  list  all  these  axioms  separately,  nor  to  cite  them  when  constructing  the  proofs  in 
section  7. 

6.1.9.  Standard  Coordinate  System. 

Definition  6.1.9.1:  The  constant  "xorigin"  denotes  a  particular  point  in  /?\ 

Definition  6.1.9.2:  The  function  ■■coordinate(X,U)"  gives  the  coordinate  of  .V  on  an  axis 
in  the  V  direction  through  the  origin. 

coordinate(X,V)  =   (X- xorigin)^ 

Definition  6.1.9.3:  The  constant  "vup"  is  the  unit  vector  in  the  upward  direction. 
length(vup)  =   1 

Definition  6.1.9.4:  The  function  "height(X)"  maps  a  point  X  into  its  height. 
height(X)  =  coordinate(X,vup) 

Theorem  6.1.9.5:  The  distance  between  rwo  points  is  greater  than  the  difference  of 
their  heights 

distance(Xl,X2)  a  height(Xl)  -  height(.V2) 

Theorem  6.1.9.6:  If  XI  is  beyond  X2  in  the  vup  direction,  then  the  height  of  XI  is 
greater  than  that  of  X2. 

beyond(Xl,X2,vup)  <=.  height(Xl)  >  height(X2) 

6.2.    Chonicles,  Motions,  and  Pseudo-Objects 

Definition  6.2.1:  The  function  'interval(Tl.r2)"  maps  two  times  7"1,  T2  onto  the  closed 
interval  from  T\  to  T2.  (72  may  be  infinity.) 

r€interval(ri,r2)  c^  T\^T<T2. 

Definition  6.2.2:  The  function  "start(rr)"  is  the  lower  bound  of  an  interval  TT .  The 
function  "end(r7")"  is  the  upper  bound  (possibly  infinity). 

start(interval(n,r2))  =  71 
end(interval(n,72))  =   72 

Definition  6.2.3:  The  function  "scenes(C)"  maps  a  chronicle  C  onto  the  set  of  scenes  in 
C. 

S  €  scenes(C)  o  3  ,-  -^  =  scene(C  .7) 
Definition  6.2.4:  The  function  "startsceDe(C)"  maps  a  chronicle  C  onto  the  first  scene  in 
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startscene(C)    =   bcene(C  ,start(duration(C  )) ) 

Definition  6.2.5:  Ihe  function  "subchronicle(C,7'r)"  maps  a  chronicle  C  and  an  time 
interval  TT  within  the  duration  of  C  into  the  portion  of  C  occuring  during  that  subchronicle. 

[rr- interval(ri.r2)  D  duration(C)  ]  => 
[  TT  =   duration(subchronicle(C,7r))  ' 
V^err  scene(C,7)  =  scene{subchronicle(C.r7'),7")  ] 

Definition  6.2.6:  The  predicate  "tail(Cl,C2)"  means  that  CI  is  the  ending  part  of  C2. 

tail(Cl,C2)  c^ 

[  3x|  rU  duration(C2)  a 

Cl  =  subchronicle(C2,intei^'al(ri.end(duration(C2))))  ] 

Definition  6.2.7:  The  predicate  "infinite(C)"  means  that  C  extends  infinitely  into  the 
future. 

Definition  6.2.8:  The  predicate  'standard-position(0,5)"  holds  if  mapping(5,0)  is  the 
identity. 

Axiom  6.2.9:  standard-positionfO  ,5')  =>  placefO,5)  =   shape((9) 

Definition    6.2.10:    The    predicate      motionlesslO.^)"    means    that    object    O    has    zero 

velocity  in  scene  S. 

motionless(0,5)  <=>  prior-velocitylS'.O)  =  0 

.^xiom  6.2.11:  If  an  object  is  motionless  during  a  chronicle,  then  it  stays  in  the  same 
place  throughout  that  chronicle. 

[  V5  5€scenes(C)  ^  motionless(0,C)  ]  => 

[  V,,  ^,  51€5cenes(C)   •  52$scenes(C)  :=.  mapping(51,0)  =  mapping(52,0)  ]  ] 

Definition  6.2.12:  The  predicate  "steady-state(C)"  means  that  everything  is  motionless 
in  C. 

steady-state(C)  «  V  >  ^[5  <;  scenes(C)  ^  motionless(0,5)  ] 

Definition  6.2.13  The  predicate  "ps-point(2 1.02)"  holds  of  a  pseudo-object  Ql  and  a 
(pseudo-)object  Q2  if  the  shape  of  Ql  is  a  single  point,  contained  in  the  shape  of  Q2,  and  Ql 
and  Q2  have  the  same  source. 

ps-point(2  l,Q2)  ^ 

[  single-point(shape(2  1))  ^  shape(Cl)  D  shape(C2)  a 
source(21)  =  source(02)  ] 

Definition  6.2.14:  The  predicate  " c\iis{Q  ,Q F ,C )"  means  that  (pseudo-)object  Q  goes 
outward  through  directed  face  QF  during  chronicle  C . 

.\xiom  6.2.15:  If  Q  exits  QF  during  C,  then  every  point  in  Q  is  weakly  in  front  of  QF 
at  some  point  during  C. 

exits(Q,QF,C)  a  ps-poiQt(QP ,Q)  => 

Bj  [56scenes(C)  a  weak-in-front(place((2/',5'),place(eF,5))  ] 

6.3.    Physics 


6.3.1.    Basic  Terms 

Definition    6.3.1.1:    The    function    "center-massfO)"    is   the   center   of   mass    (a   pseudo- 
object)  of  an  object  O. 

source-of(center-mass((?'))   =   O 

Axiom  6.3.1.2:  The  center  of  mass  of  an  object  is  inside  its  convex  hull. 
shape(center-mass(0))  €  interior(convex-hull(shape(C>))) 

Definition  6.3.1.3:  The  function  "mass(O)"  maps  an  object  O  into  its  mass. 
Axiom  6.3.1.4:  All  objects  has  positive  mass. 
mass(O)  >  0. 

Definition  6.3.1.5:  The  predicate  "fixed(O)"  means  that  object  O  is  considered  a  fixed 
obstacle.  The  opposite  is  mobiie(C>). 

mobile(O)  <=>  not(fixed((9)) 

[  fixed(O)      S€scenes(C)  ]  =>  motionlessCO.^) 

Definition  6.3.1.6:  The  predicate  "isolatedlOCC)"  means  that  no  mobile  object  in  the 
set  of  objects  00  comes  into  contact  with  any  object  outside  00  in  the  course  of  C . 

isolated(C>0,C)  «■ 
[  Vp^oo  mobile(O)  => 
[  ^o:  5.Y.V  [  5€scenes(C)  a  abut(place(0,5'),place(02,5),XX)  ]  =>  02^00  ]] 

Definition  6.3.1.7:  The  constant  "oground"  is  the  ground. 
Axiom  6.3.1.8:  The  ground  is  fixed  in  standard  position  in  anv  chronicle. 
fixed(oground)  a  standard-position(oground.5) 

Axiom  6.3.1.9:  All  objects  have  ordinary  shapes. 
ordinary(shape(0)) 

Axiom  6.3.1.10:  The  ground  is  the  only  unbounded  object. 
O^oground  =>  bounded(shape(0),Vj 

6.3.2.    Energy 

Definition  6.3.2.1:  The  function  "energy(00,5)  maps  a  set  of  objects  00  in  a  scene  S 
into  the  total  mechanical  energy  of  the  objects.  The  function  "potential-energy(00,5)"  gives 
the  potential  energy.  The  function  "kinetic-energy(0(9,5)"  gives  the  kinetic  energy. 

Axiom  6.3.2.2:  The  total  energy  is  the  potential  energy  plus  the  kinetic  energy. 
energy(00,S)  =  potential-energy(00,5)  -  kinetic-energy(C>0,5) 

Axiom  6.3.2.3:  The  total  energy  of  two  disjoints  sets  of  objects  is  equal  to  the  sum  of 
their  separate  energies 

001  n  002  =  0  =>  energyiOOWJOOl)  =  energy{C»01)  -  energy(002) 

Axiom  6.3.2.4:  The  potential  energy  of  a  single  object  is  the  height  of  its  center  of  mass 
times  its  mass  times  the  constant  g. 

potential-energy({0},5)  =  heightfplaceCcenter-massfO),^))  ■  mass(O)     g 
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Axiom  6.3.2.5:  The  constant  g  is  positive. 
g>0 

Axiom  6.3.2.6:  The  kinetic  energy  of  a  motionless  object  is  zero. 
motionless(0,5')  =>  kinetic-energydO},^ )  =  0. 

.Axiom  6.3.2.7  Kinetic  energy  is  non-negative. 
kinetic-energy(OO)  ^  0 

6.3.3.  Stable  scenes. 

Definition  6.3.3.1:  The  predicate  possibly-stablef^)"  means  that  5'  is  a  possibly  stable 
scene.  The  predicate  "necessarily-stable(5')"  means  that  S  is  necessarily  stable.  (See  section 
4.4  for  the  need  for  this  distinction.) 

Axiom  6.3.3.2:  A  necessarily  stable  scene  is  a  possibly  stable  scene. 
necessarily-stable(5)  =>  possibly-stablef5) 

Definition  6.3.3.3:  The  predicate  "necessarily-supported(0,5')"  (read  "object  O  is 
necessarily  supported  in  scene  5")  is  true  if  it  is  the  case  that,  were  every  object  except  O 
held  fixed  in  5,  then  O  would  definitely  not  move.  The  predicate  "possibly-supported(0,5')" 
is  true  if  O  would  possibly  not  move. 

.Axiom  6.3.3.4:  If  a  scene  is  necessarily  stable,  then  every  object  is  necessarily 
supported.  If  it  is  possibly  stable,  then  every  object  is  possibly  supported.  (The  converses  are 
not  true.    See  figure  26) 

necessarily-stable(5)  =>  necessarily-supported(0,5) 
possibly-stable(5)  =>  possibly-supported(C>,5) 

Axiom  6.3.3.5:  An  object  is  not  possibly  supported  unless  it  has  some  point  of  contact 
with  another  object  where  its  surface  has  a  non-positive  upward  component. 

possibly-supportedfCi')  => 

3p2..v.Y.-i,.v  [  abut(place(0,5),pIace(02,S),A'X/^)    '  X^XXA  a 

surf-norm(place(C),5),X)  ■  vup  <  0.  ] 

Axiom  6.3.3.6:  If  a  spherical  object  abuts  a  single  other  object  at  a  single  point,  then  it 
is  supported  (necessarily  or  possibly)  just  if  it  has  a  downward  surface  normal  at  the  point. 
(Note  that  this  axiom  is  the  only  one  in  our  system  which  rests  on  the  assumption  of  uniform 
mass  distribution.) 

[  sphere(shape(01))  a 

[  ^OP.xx-i  abut(place(01,5),place(0/',5),A:A:A)  «.  [OP  =  02  '  XXA  =  (XA}]]]  => 
[  [necessarily-supported(01,5)  v  possibly-supported(01,5]  o 

surf-norm(place((?/',5),Xi4)  =  -vup] 

6.3.4.  Physically  possible  chronicles 

Definition  6.3.4.1:  The  predicate  "phys-poss(C)"  means  that  chronicle  C  is  physically 
possible. 

Axiom  6.3.4.2:  Any  subchronicle  of  a  physically  possible  chronicle  is  itself  physically 
possible. 

phys-poss(Cl)  A  C2=  subchronicle(Cl, rr)  =>  phys-poss(C2) 
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Axiom  6.3.4.3:  No  two  objects  overlap  in  a  physically  possible  chronicle. 
[  phys-poss(C)   ■  5<;Scenes(C)  )  =>  not(overlap(place(01,5),place((92,S))) 

Axiom  6.3.4.4:  The  energy  of  an  isolated  system  never  increases. 

[  phys-poss(C)   •   isolated(00,C)   •   7"1<7"2  ]=> 
energytOO.scenelCT"!))  2  energy(00,scene(C,7"2)) 

Axiom  6.3.4.5:  An  infinite  chronicle  always  ends  in  a  steady  state. 

[  phys-poss(C)  and  mfinite(C)  ]  =>  B,^,  tail(Cl,C)  and  steady-state(Cl) 

Axiom  6.3.4.6:  In  a  steady-state  chronicle,  every  scene  is  possibly  stable. 

[  phys-poss(C)  and  steady-state(C)  and  5€scenes(C)  ]  => 
possibly-stable(5) 

Axiom  6.3.4.7:  If  an  object  is  at  first  inside  a  tube  and  is  later  outside  the  tube,  then  it 
must  e.xit  one  of  the  caps.t 

[  phys-poss(C)       tubetshape(Or),shape(27"0/'),shape(Q80r))  a 
place(0,51)  C  tube-inside(place('0r,51),  place(!3rO/',51),  place((2flOr,51))   \ 
disjoint  (place(0,52),  tube-inside(place(C>r,52),  place((27"0/',i'2),  place((2flOr,52))  A 
souvce(QTOP)  =  OT  and  source(efi07")  =  OT  a 
51  =  scene(C,ri)   -  52  =  scene(C,r2)  and  TKT2  ]  => 

[cxns(O.QTOP,C)   ■  sxhi,(0,QBOT ,C)  ] 

7.    Analysis  of  an  Example. 

7.1.    Problem  statement 

Consider  a  spherical  die.  and  a  radially  symmetric  funnel.  Assume  that  the  inner  radius 
of  the  funnel  is  greater  than  the  radius  of  the  die;  and  that  the  inner  side  of  a  radial  cross 
section  of  the  funnel  is  convex.  If  the  die  is  released  inside  the  funnel,  and  the  funnel  is  held 
fixed  far  from  the  ground,  then  the  die  will  eventually  fall  out  the  bottom  of  the  funnel. 

Constants  of  the  examples: 

odie  —  the  die 

ofunnel  —  the  funnel 

c  —  the  chronicle 

xx-pfunnel  —  the  planar  form  from  which  the  funnel  is  generated. 

xx-center-line  —  the  axis  of  the  funnel. 

xcenter  —  a  point  on  the  axis  of  the  funnel. 

Assumption  7.1.1:  The  die  is  a  sphere, 
sphere (shape (odie)) 

Assumption  7.1.2:  The  die  is  not  fixed. 
mobile(odie) 

Assumption  7.1.3:  The  funnel  is  the  solid  of  revolution  of  xx-pfunnel  around  xx- 
center-line. 


This  axiom  could,  and  should,  be  recast,  as  a  theorem  derived  from  (i)  a  geometrical  axiom  stating  that 
a  connnuous  curve  trom  inside  a  tube  to  outside  it  must  either  pass  through  the  tube,  through  the  top, 
or  through  the  bottom;  (li)  a  theorem  of  motion,  stating  that  any  point  in  an  objea  moving  throueh 
ame  traces  a  continuous  curve;  and  (lii)  axiom  6  3.4.?,  staung  that  objeas  do  not  overlap. 
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shape(ofunnel)  =  solid-of-revolution(xx-pfunnel,xx-center-line) 

Assumption  7.1.4:  xx-pfunnel  is  a  radial  cross  section  of  the  funnel. 
planar(xx-pfunnel  |^   xx-center-line) 

Assumption  7.1.5:  The  inner  boundary  of  xx-pfunnel  with  respect  to  xx-center-line  is 
convex. 

convex-side  (inner-dside(xx-pfunnel.xx-centerline),xx-pfunnel) 

.\ssumption  7.1.6:  The  radius  of  the  funnel  is  greater  than  the  radius  of  the  die. 
distance(xx-pfunnel,xx-centerline)  >  radius(odie)  >  0. 

Assumption  7.1.7:  The  axis  of  the  funnel  is  vertical, 
xx-centerline  =  make-line(xcenter,vup) 

Assumption  7.1.8:  The  funnel  is  oriented  in  standard  position, 
standard- position  ( of  unnel,startscene(c)) 

Assumption  7.1.9:  The  funnel  is  fixed 
fixed(ofunnel) 

Assumption  7.1.10:  The  die  is  isolated  from  everything  but  the  funnel  and  the  ground.* 
isolate  d((odie, of  unnel.oground},c) 

Assumption  7.1.11:  The  funnel  is  more  than  the  diameter  of  the  die  above  the  ground. 

XF€shape(ofunnel)  a  XG  €shape(oground)  =>  height(Xf )  -  height(XC)  >  diameter(odie) 

Assumption  7.1.12:  The  chronicle  is  eternal. 
infinite(c) 

Assumption  7.1.13:  The  chronicle  is  physically  possible. 
phys-poss{c) 

Assumption  7.1.14:  The  die  starts  from  rest. 
mot  ion  less  (odie,stanscene(c)) 

Assumption  7.1.15:  The  die  starts  from  inside  the  funnel. 

place(odie,startscene(c))  C 

tube-inside(shape(ofunnel),  s-tube-iop(shape(ofunnel)),  s-tube-bot(shape(ofunnel))) 

Prove:  The  die  exits  the  bottom  of  the  funnel. 

exits(odie,pseudo-object(ofunnel,s-tube-bot(ofunnel,vup)),c) 


'  We  need  the  ground,  because  otherwise  the  hypotheses  are  inconsistent  with  the  axioms.  The  a.xioras 
assert  that  an  infinite  chronicle  must  come  to  an  end  in  a  steady  state.  Since  we  will  show  that  there  is 
no  steady  state  for  the  die  inside  the  funnel,  we  must  provide  it  with  the  ground  to  rest  on. 
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7.2.    Analysis 

Lemma  7.2.1:  An  object  inside  a  tube  must  eventually  either  find  a  stable  state  in  the 
tube,  or  exit  one  of  the  faces. 

f  phys-poss(Cl    ■   tube(shape(07).shape(C7"OP),shape(efiOr))  A 
place(0,511  C  tube-inside(place(C'7',i'l),  place(07"O/' Jl),  place((2fl07",51))   ~ 
somcciQTOP)  =   OT  and  source(^fiOr)  =   OT  '  51  =  startscene(C)  and  infinite(C)  ] 

[  exiKiO.QTOP.C)      exns{O.QBOT,C)   . 
[  3  ,;  5€  scenes(C)   "   possibly-siable(5)   ■ 
not(disjoini(place(0,5),tube-inside(place(Or,5),place(e7"OP,5)  p\acc{QBOT ,S)))]  ] 

Proof:  Since  C  is  infinite,  it  must  end  in  some  stable  state  (6.3.4.5).  If  the  object  ever 
comes  outside,  then  it  must  exit  one  of  the  faces  (6.3.4.7).  If  it  never  comes  outside,  then  U 
is  not  outside  in  the  stable  state;  that  is,  it  is  not  disjoint  from  the  inside. 

Lemma  7.2.2:  The  energy  of  an  isolated  object  O  never  increases. 

[  isolated((90.C)   ■  5escenes(C)       phys-poss(C)  '.  "^ o\iOO  [mobile(Ol)  =>  01=0]  ]  => 
energy({0[,5 )  -^  energy({0},startscene(C)) 

Proof:  Since  the  set  00  is  isolated  in  C,  its  energy  never  increases  (6.3.4.4,  6.2.4, 
6.2.2,  6.2.1).  All  the  objects  in  00  other  than  O  are  fixed;  that  is,  they  are  always 
motionless  (6.3.1.5)  and  in  the  same  place  (6.2.11).  Their  kinetic  energy  is  thus  always  zero 
(6.3.2.6)  and  their  potential  energy  always  constant  (6.3.2.4,  6.3.1.1,  5.12),  so  their  total 
enerty  remains  constant  (6.3.2.2).  Thus  O  is  the  only  object  whose  energy  may  change. 
Since  the  energy  of  a  set  is  the  sum  of  the  individual  energies  (6.3.2.3),  the  energy  of  O 
never  increases. 

Lemma  7.2.3:  An  isolated  object  0  starting  from  rest  never  has  its  center  of  mass 
higher  than  in  the  starting  scene. 

[  isolated{OO.C)   ">  i'6scenes(C)       phys-poss(C)  /* 

'^ 0]-oo  [mobile(Ol)  =>  01  =  0]  ■  motionless(0,startscene(C))  ]  => 
height(place(center-mass(0).5))  s  height(place(center-mass(0),startscene(C))) 

Proof:  The  energy  of  O  cannot  increase  (7.2.2).  Since  it  is  motionless  in  the  starting 
scene,  its  total  energy  is  equal  to  its  potential  energy  (6.3.2.2.  6.3.2.6),  which  is  an 
increasing  function  of  the  height  of  its  center  of  mass  (6.3.2.4,  6.3.2.5,  6.3.1.5).  If  the  center 
of  mass  were  to  become  higher,  this  would  create  a  higher  potential  energy.  Since  kinetic- 
energy  IS  non-negative  (6.3.2.7),  the  total  energy  would  likewise  be  higher  (6.3.2.2),  but  this 
is  a  contradiction. 

Lemma  7.2.4:  An  isolated  convex  object  which  starts  from  rest  below  a  fixed  upward 
vertical  face  cannot  exit  through  that  face. 

[  isolated(00,C)  ■  phys-poss(C)  a  V^,.^,.,  [mobile(Ol)  =>  01  =  0)  a 

motionless(0,startscene(C))  a  convex(O)  a 

shape(QV)  =  const-face(A'A:v,vup)  a  fixed(source(2V))  a 

standard-position(source(^V),startscene(C))  a 

beyond(place(2V,startscene(C)),place(0,startscene(C)),vup)  ]  => 
not(exnsiO,QV,C)) 

Proof:  Since  0  is  convex,  the  center  of  mass  of  O  is  inside  O.  (6.3.1.2,  6.1.2.18, 
6.1.2.7).  Since  the  face  QV  is  beyond  O  in  the  starting  scene  (assumption),  every  point  in  QV 
is  higher  than  the  center  of  mass  of  O  (6.1.1.6).  If  O  were  to  exit  QV  in  C,  then  each  point 
of  O  would  have  to  be  weakly  in  front  of  QV  in  some  scene  of  C  (6.2.15).  In  particular,  the 
center  of  mass  of  O  would  have  to  be  in  front  of  QV  (6.2.13).  Since  the  source  of  QV  is 
fixed,  and  in  standard  position  facing  upwards,  (assumption),  being  weakly  in  front  of  QV 
means  being  at  least  as  high  as  QV  (6.1,515.  6.1.1.7,  6.1.1.6,  6.1.9,6).  Thus,  the  center  of 
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mass  of  O  would  have  to  be  as  high  as  than  QV ,  which  is  higher  than  the  starting  position  of 
the  center  of  mass  of  O.  (assumption)  But  this  is  impossible  by  lemma  7.2.3.  Thus,  O  does 
not  exit  QV  in  C. 

Lemma  7.2.5:  The  funnel  in  our  example  is  an  s-tube. 
s-tube{shape(ofunnel)) 

Proof:  The  funnel  is  the  surface  of  revolution  of  xx-pfunnel  around  xx-center-line 
(7.1.3).  xx-pfunnel  is  some  finite  distance  from  xx-center-line  (7.1.6);  hence  they  do  not 
overlap  (6.1.1.3,  6.1.1.4).  Hence,  the  funnel  is  an  s-tube  (6.1.7.13). 

Lemma  7.2.6:  In  our  example,  odie  does  not  exit  the  top  of  ofunnel  in  c. 
not(exits(odie,pseudo-object(ofunnel,s-tube-top(ofunnel,vup)),c) 

Proof:  We  must  show  that  the  conditions  in  (7.2.4)  hold,  with  C  =^  c.  00  =  {odie, 
ofunnel,  oground},  O  =  ofunnel,  and  QV  =  pseudo-object(ofunnel,  s-tube-top(ofunnel, 
vup)).  That  00  is  isolated  in  c  is  7.1.10.  That  c  is  physically  possible  is  7.1.13.  That  odie  is 
the  only  mobile  object  in  00  follows  from  7.1.9  and  6.3.1.8.  That  odie  is  motionless  at  the 
beginning  of  c  is  7.1.14.  That  odie  is  convex  follows  from  the  fact  that  it  is  spherical  (7.1.1, 
6.1.7.4,  6.1.2.20).  That  the  top  of  the  funnel  points  up  follows  from  7.2.5,  7.1.7,  and 
6.1.7.15.  That  the  top  of  the  funnel  has  a  fixed  source  follows  from  5.1.12.  That  the  top  of 
the  funnel  is  in  standard  position  in  the  startscene  follows  from  7.1.8  and  5.1.12.  That  the  top 
of  the  funnel  is  above  the  starting  place  of  the  die  follows  from  7.1.15,  6.1.7.15,  and 
6.1.6.11.  Thus  7.2.4  applies,  and  the  die  cannot  exit  the  top. 

Lemma  7.2.7:  The  die  does  not  ever  abut  the  top  of  the  funnel. 

5€scenes(c)  a  abut(place(odie,i'),place(ofunnel,.S),.YX/4)  => 
disjoint(XX/4  ,top(place(ofunnel,5'),vup)) 

Proof:  The  funnel  and  the  die  are  everywhere  smooth  (6.3.1.9,  6.1.5.3)  and  bounded. 
(6.3.1.10,  6.1.2.27)  Therefore,  the  surface  normal  to  the  funnel  at  a  top  point  is  vertically  up 
(6.1.5.8,  7.1.7).  Hence,  if  the  die  abuts  the  top,  the  surface  normal  to  the  die  is  vertically 
down  (6.1.5.6).  Since  the  die  is  convex,  a  downward  surface  normal  only  occurs  at  the 
bottom  of  the  die  (6.1.5.9).  Hence,  all  of  the  die,  including  the  center-of  mass,  is  higher  than 
the  contact  point  (6.1.2.26).  The  remainder  of  the  proof  is  as  in  lemma  7.2.4. 

Lemma  7.2.8:  Let  XXS  be  a  sphere.  Let  XXF  be  an  s-tube  generated  from  a  planar 
figure  with  a  convex  inner  boundary,  with  a  radius  greater  than  the  radius  of  the  sphere.  If 
XXS  abuts  the  interior  surface  of  XXF ,  it  does  so  in  a  single  point. 

[  sphere(A:X5)  a  s-tube(X.\:r)  .    A:xr=  soHd-of-revolution(A:A:f  ,XXZ.)  A 
distance(A:A:r.XXI)  >  radius(A'A:S)  a  plun^riXXF  \JXXL)  a 
convex-side(inner-dside(XXf",A'A'Z.),XXf )  a  Sibut{XXS ,XXT ,XXA)  a 
XXA  C  base(inner-boundary(A:A:r))  ] 

=>  single-point(A!X/4 ) 

Proof:  Being  a  sphere,  XXS  is  a  solid  of  revolution  around  any  line  through  its  center 
(6.1.7.11)  and  in  particular  through  the  line  through  its  center  parallel  to  XXL  (6.1.2.4).  Call 
this  line  XXLl.  If  these  two  lines  were  the  same,  then  XXS  could  not  abut  XXT  (6.1.2.12), 
since  all  points  of  XXS  are  within  radius(XX5)  of  its  center  (6.1.7.6),  and  hence  within 
radius(XX5)  of  XXLl  (6.1.1.3),  while  all  points  of  XXT  are  further  than  radius(XX5)  from 
XXL  (assumption).  Therefore,  XXL* XXLl  so  XXA  lies  in  a  plane  XXP  which  includes  XXL 
and  XXLl  (6.1.7.12).  XXA  therefore  lies  m  the  intersection  of  XXT  and  XXP,  and  in  the 
intersection  of  XXS  with  XXP .  The  intersection  of  XXT  and  XXP  consists  of  two  cross-sections 
which  are  separated  by  the  inner  diameter  of  the  funnel  (7.2.5,  6.1.7.16).  Since  the  inner 
diameter  of  the  funnel  is  greater  than  the  diameter  of  the  sphere,  the  sphere  can  intersect 
only   one   of   these   cross   sections.    (7.1.6,   6.1.7.7,   6.1.2.23)    The    intersection   of   this  cross 
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section  with  the  inner  boundary  of  the  funnel  is  congruent  to  the  inner  side  of  XXF 
(6.1.7.17),  and  is  therefore  a  convex  curve  (7.1.5).  Since  the  intersection  of  the  sphere  with 
the  plane  contains  the  center,  which  is  not  a  boundary  point  (6.1.7.9),  this  intersection  does 
not  consist  of  a  single  boundary  point,  and  is  therefore  a  circle  (6.1.7.8),  which  is  strongly 
convex  (6.1.7.2).  Since  these  two  planar  cross  sections  abut  (6.1.3.8)  they  must  abut  in  a 
single  point  (6.1.4.6),  giving  the  desired  result. 

Lemma  7.2.9:  If  a  convex  shape  abuts  the  inside  of  a  p-tube,  then  any  abutment  points 
of  the  shape  with  the  p-tube  lies  on  the  top,  or  the  bottom,  of  the  p-tube. 

[  p-tubt(XXT,P\,P2,V)  A  con\ex(XXO)  a 

abutiXXO ,tube-mside(XXT ,P  1.P2))  a  abux(XXO,XXT ,XXA)  ]  => 
[  XXA  D  top(A:A'7",V)  V  XXA  D  bottom (XXT",!/)  ] 

Proof:  Since  XXO  does  not  overlap  XXT  (6.1.2.14),  it  cannot  abut  the  inside  of  XXT  at 
an  interior  point  of  the  common  boundary  of  XXT  with  its  inside  (6.1.5.7,  6.1.5.5,  6.3.1.9, 
6.1.5.3),  which  is  the  inner  surface  (6.1.6.8).  Therefore  XXO  must  abut  the  inside  of  XXT  on 
either  PI  or  P2  (6.1.6.7,  6.1.6.8,  6.1.6.9).  Hence,  XXO  must  be  weakly  in  front  of  either  P\ 
or  P2  (6.1.5.16,  6.1.5.17).  But  the  only  points  in  XXT  which  are  weakly  in  front  of  PI  (P2) 
are  the  top  (bottom)  of  XXT  (6.1.6.10,  6.1.3.10,  6.1.2.26).  Hence,  the  abutment  points  must 
lie  in  the  top  (bottom)  oi  XXT . 

Lemma  7.2.10:  If  any  part  of  the  object  odie  is  in  the  inside  of  the  funnel  ofunnel,  then 
odie  is  either  disjoint  from  ofunnel  or  odie  abuts  ofunnel  in  a  single  point  in  the  internal 
surface. 

[  5€scenes(c)  a  not(disjoint(place(odie,5),place(s-tube-inside(ofunnel),5)))  ]  => 
[  [3x  abut(place(odie,5),place(ofunnel,5'),(X})  ,'  X€  place(inner-boundary(ofunnel),5)  ]  v 
disjoint(place(ofunnel.5'),place(odie,5' ))  ] 

Geometrically,  odie  and  ofunnel  must  either  be  disjoint,  abut,  or  overlap  (6.1.2.15, 
6.3.1.9,  6.1.5.3);  physically,  they  cannot  overlap  (6.3.4.3,  7.1.13).  Since  odie  is  convex 
(7.1.1,  6.1.7.4,  6.1.2.20),  and  ofunnel  is  a  p-tube  (7.2.5,  6.1.7.15),  if  any  part  of  odie  is  in 
the  inside  of  ofunnel,  then  all  abutment  points  of  odie  with  ofunnel  are  in  the  inner  boundary 
of  ofunnel  (7.2.9).  But  odie  and  ofunnel  satisfy  the  conditions  in  (7.2.8),  by  virtue  of  7.1.1, 
7.2.5,  7.1.4,  7.1.5,  and  7.1.6.  Hence,  there  is  only  one  abutment  point  between  odie  and  the 
inner  boundary  of  ofunnel,  which  proves  the  result. 

Lemma  7.2. II:  If  odie  abuts  the  ground,  then  it  is  disjoint  from  the  inside  of  the 
funnel. 

abut(odie,oground,XXA  )  =>  disjoint(odie,s-tube-inside(ofunnel)) 

Proof  If  odie  abuts  the  ground,  then  any  point  in  odie  is  within  diameter(odie)  of  the 
ground  (6.1.2.23,  6.1.2.12).  Therefore,  the  height  of  any  point  in  odie  minus  the  height  of 
any  point  in  the  ground  is  at  most  the  diameter  of  odie  (6.1.9.5).  But  all  points  in  ofunnel  are 
more  than  diameter(odie)  above  any  point  in  the  ground  (7.1.11),  and  all  points  in  the  inside 
of  ofunnel  are  above  the  bottom  of  ofunnel  (6.1.6.11).  Hence,  the  inside  of  ofunnel  must  be 
disjoint  from  odie. 

Lemma  7.2.12:  No  point  on  the  inner  surface  of  a  smooth  s-tube  constructed  from  a 
planar  figure  with  a  convex  inner  boundary  points  vertically  upward. 

[  s-tube(XXr)  A  A'X7=solid-of-revolution(A:XF,XA:L)  a 
planzT(XXF\jXXL)  a  convex-side(inner-side(X>rf,XXZ.),XXF)  a 
smoothiXXT J()  a  X€inner-boundary(XXr)  ]  => 

surf-norm(XA:r,X)  i=  vup 

Proof:  Since  XXT  is  a  solid  of  revolution,  the  surface  normal  to  XXT  at  X  lies  in  the 
plane  of  X  and  XXL  (6.1.7.18),  and  is  equal  to  the  normal  to  the  cross  section  of  XXT"  by  that 
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plane  which  contains  X  (6.1.5.10).  This  cross-section  is  congruent  to  XXF  under  a  mapping 
that  leave  XXL  fixed  (6.1.7.16);  by  the  invariance  of  geometric  properties  under  mappings, 
its  inner  boundary  is  likewise  convex,  and  so  has  no  normal  in  the  direction  of  XXL,  which  is 
vup  (6.1.4.8). 

Lemma  7.2.13:  A  convex  mobile  object  which  abuts  only  the  bottoms  of  other  objects  is 
not  possibly  supported. 

[  convex(O)  '  mobile(O)  a 

^  OP  v.v-i  [abut(place(0,5),place(0/',5),XX/\)  =>  XXA  C  bottom(0/',vup)  ]  ]  => 
no  t(  possibly -supported(  0,5)) 

Proof:  The  surface  normal  to  any  of  these  other  objects  OP  at  their  bottom  is  vertically 
down  (6.3.1.9,  6.1.5.3,  6.1.5.8);  hence,  the  surface  normals  to  O  there  are  all  vertically  up 
(6.1.5.5,  6.1.5.6);  hence,  they  do  not  have  negative  dot  products  with  vup  so  that  O  is 
unsupported  (6.3.3.5). 

Lemma  7.2.14:  In  any  stable  scene,  odie  is  disjoint  from  the  inside  of  ofunnel. 

[  possibly-stable(5)  a  5€scenes(c)  ]  => 
disjoint(place(odie,5),  place  (s- tube-inside  (of  unnel),5)) 

Proof:  If  odie  is  not  disjoint  from  the  inside  of  ofunnel,  then  it  does  not  abut  the  ground 
(7.2.11);  since  it  is  isolated  from  everything  but  the  ground  and  the  funnel,  it  does  not  abut 
anything  but  the  funnel  (7.1.10,  7.1.2,  6.3.1.6).  Furthermore,  if  it  is  not  disjoint  from  the 
inside  then  any  abutment  points  are  either  on  the  inner  boundary,  or  on  the  top,  or  on  the 
bottom  (7.2.9,  7.2.10).  It  cannot  attain  a  position  where  it  abuts  the  funnel  on  top  (7.2.7).  If 
it  abuts  on  the  inner  boundary,  then  it  is  abutting  in  a  single  point  (7.2.10),  with  a  non- 
vertical  surface  normal  (7.2.12);  hence,  the  surface  normal  of  the  die  cannot  be  vertically 
downward  (6.1.5.5,  6.1.5.6).  The  die,  being  spherical  (7.1.1).  can  only  be  supported  at  a 
single  point  if  its  surface  normal  is  downward  at  that  point  (6.3.3.6).  Thus,  it  is  not 
supported  on  the  inner  boundary.  If  it  abuts  on  the  bottom,  then  it  is  unsupported,  by 
(7.2.13).  If  it  does  not  abut  at  all,  then  it  is  certainly  unsupported  (6.3.3.5).  Since  it  is 
unsupported,  the  scene  is  unstable  (6.3.3.4). 

Theorem  7.2.15:  The  die  comes  out  the  bottom  of  the  funnel. 

exits(odie,s-tube-bottom(ofunnel),c) 

Proof:  Since  c  is  infinite  (7.1.12),  it  must  end  in  some  possibly  stable  state  (7.1.13, 
6.3.4.5.  6.3.4.6,  6.3.4.2,  6.2.6).  No  stable  state  can  occur  until  the  die  is  disjoint  from  the 
inside  of  the  funnel.  (7.2  14).  Therefore,  the  die  must  either  exit  the  top  or  the  bottom 
(7.2.1).  Since  it  cannot  exit  the  top  (7.2.6)  it  must  exit  the  bottom. 

8.    Conclusions 

The  strengths  and  limitations  of  this  theory  are  both  quite  evident.  On  the  positive  side: 
We  have  given  a  formal  analysis  of  a  class  of  problems  far  beyond  the  scope  of  any  previous 
AI  theory,  and  we  have  completed  the  analysis  within  the  boundaries  of  first-order 
deduction.  Our  analysis  suggests  that  a  qualitative  physics  for  solid  objects  should  include  the 
following  features,  among  others: 

•  A  rich  geometrical  theory,  including  topological,  metric,  and  differential  descriptors, 
and  special  shapes.  (The  concepts  of  point-set  topology  used  in  section  6,  such  as 
interior  and  boundary,  would  probably  not  require  explicit  representation  in  a 
computational  system.) 

•  An  account  of  the  behavior  of  physical  systems  over  extended  intervals  of  time.  Such  an 
account  should  incorporate  constraints  placed  by  one  object  on  another;  conservation 
laws,  especially  conservation  of  energy;  the  principle  that  a  physical  system  tends 
towards   a  stable   resting   point;  and   an  account  of  the   net  effects  of  collisions  over 
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extended  time  periods.  As  regards  collisions,  the  complete  mathematical  and  physical 
analysis  is  not  obvious.  It  may  be  possible  to  carry  out  such  a  complete  analysis;  or  it 
may  be  possible  to  rest  on  heuristic  rules  such  as  suggested  in  section  3. 

•  The  ability  to  determine  the  existence  of  a  stable  configuration  of  objects  within 
qualitatively  described  geometrical  constraints. 

•  The  ability  to  calculate,  exactly  or  qualitatively,  important  physical  parameters  such  as 
the  center  of  mass.  (Bundy  and  Byrd  (1983)  present  some  interesting  techniques  for  this 
kind  of  problem.) 

•  The  ability  to  calculate  bounds  on  the  effect  of  small  perturbations  on  physical 
parameters  or  physical  behavior. 

The  analyses  in  section  3  use  all  these  features.  The  language  and  axioms  in  section  6 
includes  examples  of  them  all  except  the  analysis  of  collisions  and  of  perturbations. 

On  the  negative  side;  We  have  not  shown  that  this  type  of  analysis  is  extensible  to  cover 
all,  or  most,  qualitative  reasoning  in  this  domain.  We  have  not  shown  that  such  an  extension 
would  be,  in  the  long  run,  any  more  parsimonious  than  simply  enumerating  special  cases,  as 
in  the  rule-based  method  rejected  in  section  2.  We  have  not  shown  that  any  effective 
computational  methods  can  be  developed  on  the  basis  of  this  theory. 

Actually,  none  of  these  doubts  can  be  effectively  resolved  without  a  robust  running 
program.  The  question  of  extensbility  cannot  be  answered  without  analysing  a  fairly  large 
number  of  examples.  The  question  of  parsimony  cannot  be  resolved  without  analysing  a  very 
large  number  of  examples.  But  no  one  in  the  world  is  going  to  have  the  patience  to  analyse  a 
very  large  number  of  examples  by  hand  in  the  style  of  sections  6  and  7.  The  only  practical 
way  to  show  that  our  theory  is  robust,  even  in  a  logical  sense,  is  to  create  a  running  program, 
and  show  that  it  works  on  very  many  examples. 

The  complexity  of  our  theory  and  its  'eliance  on  formal  physics  and  geometry  will 
almost  certainly  elicit  the  following  criticism:  "You  have  shown  that  you  can  solve  common 
sense  physical  problems  by  throwing  all  the  physics  and  geometry  you  know  at  them.  So 
what?  Naturally,  classical  physics  and  mathematics  are  sufficient  to  prove  that  a  die  falls 
through  a  funnel.  Who  would  have  doubted  it?" 

There  are  two  answers  to  this.  First,  we  have  seen  that  the  analysis  is  not  simple.  For 
several  cases  in  section  3,  I  had  to  employ  high-level  principles  that  I  have  not  been  able  to 
prove  correct.  No  doubt,  this  in  part  reflects  my  ignorance  or  incompetence.  Still,  there  is 
good  reason  to  believe  that  many  physical  problems  whose  solutions  are  clear  to  a  common 
sense  understanding  are  quite  difficult  to  analyse  in  this  way.  In  view  of  this  difficulty,  the 
statement  that  the  analysis  is  even  possible  becomes  not  wholly  trivial. 

Second,  and  more  importantly,  the  details  of  the  analysis  itself  is  suggestive.  It  is 
particularly  interesting  to  see  what  kinds  of  conditions  must  be  placed  on  the  problem  for  the 
solution  to  follow  soundly.  For  example,  the  features  enumerated  above  were  extracted  by  a 
consideration  of  requirement  of  the  analysis.  For  another  example,  the  observation  that 
physical  behavior  is  sensitive  to  local  changes  in  the  surface  normal  came  as  a  surprise  to  me, 
and  I  only  realized  it  when  I  began  the  formal  analysis  of  the  problem.  This  "logical 
framework"  can  be  viewed  as  a  collection  of  useful  lemmas  from  Newtonian  mechanics,  some 
of  which  are  unproven.  It  is  of  some  interest  to  know  what  lemmas  can  be  used  for  what 
problems. 

It  is  also  important  to  note  what  concepts  have  been  unnecessary  in  solving  these 
problems.  For  the  examples  in  section  3,  we  did  not  have  to  reason  about  velocities,  except 
the  zero  velocity,  momentum,  or  angular  momentum.  If  we  may  use  the  cited  rules  which 
determine  possible  stability  directly  from  geometrical  specifications,  we  may  also  omit  any 
reasoning  about  forces  or  torques.  Thrs  is  not  to  suggest  that  these  concepts  can  be  ignored  in 
a  larger  space  of  examples,  nor  that  these  examples  could  not  be  given  an  alternative  analysis 
which  did  use  these  concepts.  Still,  it  is  important  to  know  how  powerful  a  limited  set  of 
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concepts  is  in  solving  particular  classes  of  problems. 

There  are  two  areas  in  which  I  feel  that  further  formal  analysis  would  be  valuable.  The 
first,  as  mentioned,  is  the  analysis  of  collisions.  The  second,  which  is  not  at  all  addressed  m 
this  paper,  is  the  analysis  of  problems  involving  metric  time  and  velocity.  To  deduce  that  the 
die  comes  out  the  bottom  of  the  funnel  is  not,  after  all,  all  that  useful.  We  would  like  to 
know  roughly  when  will  come  out;  next  year,  in  an  hour,  in  a  minute,  in  a  rwinkling  of  an 
eye?  We  will  also  like  to  know  something  about  its  velocity  on  exit;  to  be  sure,  for  example, 
that  it  will  not  come  hurtling  out  in  a  nearly  horizontal  direction.  Nothing  in  the  theory 
developed  so  far  addresses  either  of  these  issues. 

Our  short  term  goals  are  to  complete  these  aspects  of  formal  analysis  and  to  begin 
implementation  by  developing  an  adequate  geometric  representation  and  inference  system. 
Ultimately,  we  want  to  implement  a  physical  reasoning  system  with  all  the  feature  mentioned 
above. 
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Fiqure  1:  A  die  is  released  in  a  funnel 
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Figure  2:  DiscretizecJ  space  and  velocity  in  FROB 
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Figure  3:  One  die  inside  another  released  in  a  tunnel 


Figure  4:  A  spherical  die  inside  a  radially  symmetric  funnel 
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The  die  cannot  come  out  the  hole  because  of  the  fixed  barrier 

Figure  5 
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Torus  T  is  cut  by  plane  P 


Cross  section  of  T  hy   P 
C  is  a  curve  in  the  cross  section, 
p  is  a  point  inside  C 


Figure  6 
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Figure  /:  Boomerang  v/ith  center  of  mass 


Figure  8 :  -  Two  dice  in  a  funnel 
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The  dent  will  spin  the  ball  around  without  energy  loss. 

Figure  9 
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Little   energy    is    lost    in   a    side-sv/ipe. 
Figure    10 
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The  ball  rolls  up  .and  down  without  energy  loss. 

Figure  11 
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f  is  the  maximum  angle  between  the  normal  and  tlie  incoming  velocity 

Figure  12:  Head  on  collision 
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N  is  the  normal  to  the  surface, 
c  is  the  center  of  mass. 
1  is  a  line  through  c. 

is  the  maximum  angle  between  K  and  1 

Figure  14 :  Distorted  sphere 


Figure  15:  Spherical  die  contains  another  die. 
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Side  view: 
Q  is  the  minimuin  slope 
at  a  contact  point. 


Top  view: 
A  and  B  are  contact  points. 
Nl  and  U2   are  horizontal 
projections  of  surface  normals 
at  A  and  B.  u  is  a  vector 
between  Nl  and  N2.  ^  is  an 
upper  bound  on  the  angle  between 
u  and  Nl,  or  u  and  N2. 


Figure  16:  General  conditions  for  support 
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One  object  may  contact  another  in  a  face,  a  line,  or  a  point 

Figure  17 
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The  distribution  of  the  weight  of  .X  between  A,  B,  and  C  is  underdetermined 

Figure  18 
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The  horizontal  normal  forces  betv;een  A  anr]  B  are  underdetermined. 
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Startinq  situation  One  possible  result       Another  result 

Figure  20:  Many  possibilities  in  a  multi-object  collision 
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starting  situation 
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End  result 
if  A  hits  B  before  C. 


End  result 
if  C  hits  B  before  A. 


The  end  result  of  a  collision  may  be  discontinuous  under 
small  changes  in  timing. 

Figure  21 


The  rebound  direction  of  A  depends  on  tlie  distribution  of  force 

along  the  surface  of  contact. 


Figure  22 
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If  A  and  B  have  different  velocities,  tl.ere  is  no  consistent  solution 

Fiqnie  2  3 
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The  inner  side  of  XXP  v;ith  res{5ect  to  XXL 
is  the  part  of  its  boundary  bet-.veen  a  an(]  b. 

Fiqui-e  24 
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XXT  is  the  outer  shell. 
PI,  P2  are  the  upper  and  lov;er  faces. 
XXCl  and  XXC2  are  the  top  and  bottom 
thirds  of  the  internal  cylinder. 

Figure  25:  Definition  of  a  tube. 
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B  and  A  are  both  "possibly  supported" 
(if  either  was  fixed,  the  otlier  would  be  supported.) 
However,  if  they  are  both  mobile,  tliey  will  fall  together. 

Figure  26 
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